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ABSTRACT 

We report on the X-ray evolution over the last ^9 Gyr of cosmic history (i.e., since z = 1 .4) of late-type galaxy 
populations in the Chandra Deep Field-North and Extended Chandra Deep Field-South (CDF-N and E-CDF-S, 
respectively; jointly CDFs) survey fields. Our late-type galaxy sample consists of 2568 galaxies, which were 
identified using rest-frame optical colors and HST morphologies. We utilized X-ray stacking analyses to in- 
vestigate the X-ray emission from these galaxies, emphasizing the contributions from normal galaxies that are 
not dominated by active galactic nuclei (AGNs). Over this redshift range, we find significant increases (factors 
of ^5-10) in the X-ray-to-optical mean luminosity ratio (L x /L fl ) and the X-ray-to-stellar-mass mean ratio 
(L X /A/*) for galaxy populations selected by L B and M*, respectively. When analyzing galaxy samples selected 
via SFR, we find that the mean X-ray-to-SFR ratio (L X /SFR) is consistent with being constant over the entire 
redshift range for galaxies with SFR =1-100 M© yr _1 , thus demonstrating that X-ray emission can be used as 
a robust indicator of star-formation activity out to z ~ 1.4. We find that the star- formation activity (as traced by 
X-ray luminosity) per unit stellar mass in a given redshift bin increases with decreasing stellar mass over the 
redshift range z = 0.2-1, which is consistent with previous studies of how star- formation activity depends on 
stellar mass. Finally, we extend our X-ray analyses to Lyman break galaxies at z ~ 3 and estimate that Lx/L B 
at z ~ 3 is similar to its value at z = 1 .4. 

Subject headings : cosmology: observations — surveys — galaxies: normal — galaxies: spiral — galaxies: 
star formation galaxies: active — galaxies: — X-rays: galaxies — X-rays: general 


1. introduction 

Investigations focusing on global changes in star-formation 
activity and stellar-mass build-up in field galaxies have pro- 
vided significant insight into the physical evolution of galax- 
ies and their constituent stellar populations. It has now been 
well-established that the global star-formation rate density has 
declined by roughly an order of magnitude since z « 1-1.5 
(e.g., Lilly etal. 1996; Madau etal. 1996; Steidel etal. 1999; 
Hopkins etal. 2004; Perez-Gonzalez etal. 2005; Schimi- 
novich etal. 2005; Colbert etal. 2006). Recent investiga- 
tions into the details of this evolution have shown that the 
star- formation history of a given galaxy depends strongly on 
its stellar mass (e.g., Cowie etal. 1996; Juneau etal. 2005; 
Bundy etal. 2006; Noeske etal. 2007a, b; Zheng etal. 2007); 
the peak star-formation epoch for the most massive galaxies 
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occured earlier in cosmic history than it did for galaxies with 
lower masses. 

X-ray studies of normal late-type galaxies (i.e., those 
that are not dominated by luminous active galactic nuclei 
[AGNs]) have shown that X-ray emission provides a useful, 
relatively-unobscured measure of star-formation activity (e.g., 
Bauer etal. 2002a; Cohen 2003; Ranalli etal. 2003; Colbert 
etal. 2004; Grimm etal. 2003; Gilfanov etal. 2004a; Per- 
sic etal. 2004; Persic & Rephaeli 2007; however, see Barger 
etal. 2007). In normal galaxies, X-ray emission originates 
from X-ray binaries, supemovae, supernova remnants, hot 
(^0.2-1 keV) interstellar gas, and O-stars (see, e.g., Fab- 
biano 1989, 2006 for reviews). Sensitive Chandra and XMM- 
Newton studies of local late-type galaxies have shown that 
high-mass X-ray binaries (HMXBs) and low-mass X-ray bi- 
naries (LMXBs) typically dominate the total non-nuclear 
X-ray power output (e.g., Zezas et al. 2002; Bauer et al. 2003; 
Soria & Wu 2003; Swartz et al. 2003; Jenkins etal. 2005; Kil- 
gard etal. 2005; however, see, e.g., Doane etal. 2004). Ob- 
servations indicate that the integrated X-ray emission from 
HMXB and LMXB populations trace galaxy star-formation 
rate (SFR) and stellar mass (M*), respectively. For example, 
using Chandra observations of 32 local galaxies, Colbert et al. 
(2004) found that the summed 0.3-8 keV non-nuclear point- 
source emission from a given galaxy (L XP ) can be approxi- 
mated as L X p + SFR, where a and j3 are constants. 

Therefore, galaxies having relatively high star-formation rates 
per unit mass (specific star-formation rates, SSFRs) gener- 
ally have dominant X-ray point-source contributions from 
HMXBs (e.g., late-type star- forming galaxies), while those 
with relatively low SSFRs have point-source emission primar- 
ily from LMXBs (e.g., massive early-type galaxies). 

If the X-ray binary populations are similarly dominating the 
normal-galaxy X-ray power output over a significant fraction 


of cosmic time, then there should be a rapid increase in the 
globally-averaged X-ray luminosity of normal star-forming 
galaxies with cosmic look-back time in response to the in- 
creasing global star- formation rate density (e.g., Ghosh & 
White 2001). In this scenario, HMXBs trace the immedi- 
ate star-formation rate of a galaxy and LMXBs trace its star- 
formation history with a lag of a few Gyr. With the advent of 
deep Chandra and XMM-Newton surveys (see, e.g., Brandt 
& Hasinger 2005 for a review), it has become possible to 
study the X-ray properties of normal galaxies out to z ^ 1 
and z«0.3, respectively (see, e.g., Homschemeier et al. 2000, 
2003; Alexander etal. 2002; Georgakakis etal. 2007; Geor- 
gantopoulos etal. 2005; Lehmer etal. 2006, 2007; Kim etal. 
2006; Tzanavaris etal. 2006; Rosa-Gonzalez etal. 2007). Ev- 
idence for a global increase in the X-ray activity with red- 
shift for normal late-type galaxies has since been mounting. 
Studies of the X-ray luminosity functions of X-ray-detected 
normal galaxies have found that evolves as (1+z) 1,5-3 
over the redshift range z ~ 0-1.4 (e.g., Norman etal. 2004; 
Georgakakis etal. 2006; Ptak etal. 2007; Tzanavaris & Geor- 
gantopolous 2008). Additionally, X-ray stacking analyses 
have enabled investigations of more representative optically- 
selected galaxy populations over the majority of cosmic his- 
tory (z « 0.1-4; e.g., Brandt et al. 2001 ; Homschemeier et al. 
2002; Nandra etal. 2002; Georgakakis etal. 2003; Reddy & 
Steidel 2004; Laird etal. 2005, 2006; Lehmer etal. 2005a). 
These studies have found that the average X-ray luminosi- 
ties of normal late-type galaxies increases with redshift out 
to z ~ 1 .4-3. For example, using awl Ms exposure of a sub- 
region within the Chandra Deep Field-North, Homschemeier 
etal. (2002, hereafter H02) tentatively observed a factor of 
£s 2-3 increase in Lx /L b from z = 0 to 1.4 for L* B galaxies. 
Despite these promising initial constraints, the details of the 
X-ray evolution of normal late-type galaxy populations, in- 
cluding dependences on the physical properties of galaxies 
(e.g., optical luminosity, stellar mass, environment, and star- 
formation rate), have remained unexplored. 

In this paper, we aim to improve significantly upon con- 
straints for the X-ray evolution of normal late-type galax- 
ies (e.g., H02). We study for the first time how the X-ray 
properties of late-type field galaxies evolve as a function 
of optical luminosity, stellar mass, and star-formation rate 
over the redshift range of z = 0-1.4. We construct late- type 
galaxy samples located in two of the most well-studied ex- 
tragalactic X-ray survey fields, the Ms Chandra Deep 
Field-North (CDF-N; Alexander et al. 2003) and the Extended 
Chandra Deep Field-South (E-CDF-S), which is composed 
of the central Ms Chandra Deep Field-South (CDF-S; 
Giacconi etal. 2002) and four flanking «250 ks Chandra 
observations (Lehmer etal. 2005b). These Chandra Deep 
Fields (hereafter CDFs) reach 0.5-2 keV detection limits of 
^2.5 x 10 -17 ergs cm -2 s _1 in the most sensitive regions and 
^3 x 10~ 16 ergs cm -2 s _I over the majority of the CDFs; 
these levels are sufficient to detect moderately-powerful X-ray 
sources (Lo.5-2 keV ~ 10 41 5 ergs s _1 ) at z = 1 .4 and z = 0.6, re- 
spectively. Therefore, the CDFs comprise an unprecedented 
data set for effectively studying the X-ray emission and evolu- 
tion of cosmologically distant normal galaxies with minimal 
contamination from powerful AGNs. 

The Galactic column densities are 1.3 x 10 20 cm -2 for 
the CDF-N (Lockman 2004) and 8.8 x 10 19 cm -2 for the 
E-CDF-S (Stark etal. 1992). All of the X-ray fluxes and lu- 
minosities quoted throughout this paper have been corrected 
for Galactic absorption. Unless stated otherwise, we quote 


optical magnitudes based upon the AB magnitude system 
for HST photometry and the Vega magnitude system every- 
where else. In the X-ray band, we make use of three standard 
bandpasses: 0.5-2 keV (soft band [SB]), 2-8 keV (hard band 
[HB]), 0.5-8 keV (full band [FB]). Throughout this paper, we 
make estimates of stellar mass and star-formation rates us- 
ing a Kroupa (2001) initial mass function (IMF); when mak- 
ing comparisons between these estimates and those quoted in 
other studies, we have adjusted all values to correspond to our 
adopted IMF. H 0 = 70 km s -1 Mpc -1 , Qm = 0.3 , and f2A = 0.7 
are adopted throughout this paper (e.g., Spergel etal. 2003), 
which imply a look-back time of 7.7 Gyr at z = 1 . 

2. LATE-TYPE GALAXY SAMPLE SELECTION 

We constructed an optically selected sample of late-type 
galaxies within the CDFs to use for our subsequent anal- 
yses. We restricted our galaxy selection to regions of the 
CDFs where Hubble Space Telescope (HST) observations 
were available to allow the best possible morphological clas- 
sifications. The HST observations in the CDFs have been 
carried out via the Great Observatories Origins Deep Sur- 
vey (GOODS; Giavalisco etal. 2004a) and Galaxy Evolution 
from Morphology and SEDs (GEMS; Rix etal. 2004; Cald- 
well etal. 2005) programs; these surveys cover &90% of the 
Chandra - observed regions of the CDFs with the Advanced 
Camera for Surveys (ACS). In the GOODS and GEMS re- 
gions, photometry was available for four (£435 , V6O6, hi 5 , and 
Z850) and two (V6O6 and zsso) ACS passbands, respectively. 

2.1. Galaxy Selection Footprint 

We began building our sample by selecting all galaxies hav- 
ing zsso < 23. This initial selection criterion was motivated by 
(1) the availability of deep zsso band imaging over all of the 
CDFs, (2) the fact that the zsso emission probes mass-tracing 
rest-frame optical light redward of the 4000 A break for galax- 
ies at z £ 1 , which constitutes the large majority of our sam- 
ple, (3) the availability of reliable spectroscopic and photo- 
metric redshifts (see details below), and (4) the high multi- 
wavelength detection fractions for these sources, which allow 
us to determine informative rest-frame spectral energy distri- 
butions (SEDs). In order to isolate most effectively distant 
X-ray-detected AGNs, we further restricted our sample to in- 
clude only galaxies that were located in the most sensitive 
areas of the CDFs where the Chandra point-spread function 
(PSF) was small. We therefore chose to include sources hav- 
ing optical positions that were within 6'0 of at least one of 
the six Chandra aimpoints in the CDFs 12 ; the corresponding 
total areal footprint is «0.18 deg 2 . Futhermore, we removed 
obvious Galactic stars that were identified via optical spectral 
properties or the presence of obvious diffraction spikes in the 
Zsso band images. Under these restrictions, we found 6905 
galaxies. 

2.2. Redshifts 

We cross-correlated our initial sample of 6905 galaxies 
with available spectroscopic and photometric redshift cata- 
logs (e.g., Barger etal. 2003; Le Fevre etal. 2004; Szokoly 
etal. 2004; Wirth etal. 2004; Wolf etal. 2004; Mobasher 
etal. 2004; Mignoli etal. 2005; Vanzella etal. 2005, 2006; 
Grazian etal. 2006; Ravikumar etal. 2007; Silverman etal. 

12 For the CDF aimpoints, see Tables 1 and Al of Alexander etal. (2003) 
for the Ms CDF-N and «1 Ms CDF-S, respectively, as well as Table 1 of 

Lehmer et al. (2005b) for the ^250 ks E-CDF-S. 



Fig. 1 Rest-frame U-V color versus absolute V-band magnitude My (color-magnitude diagrams) for the 3259 z ~ 0-1 .4 galaxies with zgso < 23 that were 
within 6f0 of at least one CDF aimpoint. Each panel shows the color-magnitude relation for a given redshift bin {annotated in the upper left-hand corners ). 
The dashed line in each panel represents the estimated division between red (large U-V values) and blue (small U-V values) galaxy populations, which was 
estimated using equation 1 and the median redshift of the galaxies in each bin. Open symbols highlight red-sequence and blue-cloud galaxies that were visually 
reclassified as late-type ( open circles ) and early-type {open squares ) galaxies, respectively. The inset histogram in each panel shows the distribution of rest-frame 
U-V . The vertical dotted line in each inset plot indicates the estimated division between red and blue galaxy populations, which was calculated using equation 1 
and the median redshift and My for galaxies in each redshift bin (see § 2.3). 


2008). All galaxies that did not have spectroscopic redshifts 
were located in the E-CDF-S where highly accurate (median 
Sz/ 1 + z ~ 0.02 for galaxies with zsso < 23) photometric red- 
shifts were available via COMBO- 17 (Classifying Objects by 
Medium-Band Observations in 17 Filters; Wolf etal. 2004). 
In total, 6683 («97%) of our sources had either spectroscopic 
or photometric redshifts. Visual inspection of the 222 galaxies 
in the E-CDF-S without redshifts indicated that these sources 
were mainly faint galaxies near bright stars, as well as a hand- 
ful of sources that were subgalactic features within relatively 
nearby galaxies. 

Whenever possible, we adopted spectroscopic redshifts as 
the most accurate redshifts for our galaxies. Using the red- 
shift information, we filtered our sample to include only 
sources with z < 1 .4 in the «2 Ms CDF-N, z < 1 in the 
~1 Ms CDF-S, and z < 0.6 in the ^250 ks regions of the 
E-CDF-S; these redshift limits represent the largest distances 
at which we would expect to identify moderately luminous 


(£ 0 . 5-2 kev £ 10 41 5 ergs s' 1 ) AGNs effectively in each respec- 
tive field. In total, 3259 galaxies remained after filtering our 
sample based on redshift properties. We used spectroscopic 
redshifts for 1351 (^41%) galaxies and photometric redshifts 
for the remaining 1908 (^59%) galaxies. 

2.3. Re st- Frame Color and Morphological Selection 

The optical-color distribution for field galaxies has been 
shown to be bimodal, separating “red” and “blue” galaxy pop- 
ulations (e.g., Strateva etal. 2001; Hogg etal. 2002a; Blan- 
ton etal. 2003; Baldry etal. 2004). Studies of the color- 
magnitude relation for distant galaxy populations have shown 
that this color bimodality is observed to persist out to at least 
z ~ 1-1 .5 (see, e.g., Bell etal. 2004a; Faber etal. 2007; Labbe 
et al. 2007), thus providing an excellent quantifiable means for 
separating late-type and early-type galaxy populations. We 
therefore filtered our galaxy sample to include only sources 
that had blue rest-frame optical colors, as expected for late- 



type galaxies that contain young stellar populations. To es- 
timate rest-frame U-V colors and absolute £/, B, and V 
band magnitudes (My, M B , and My , respectively) for each 
of the 3259 galaxies in our sample, we utilized the available 
optic al/near-IR data. 

For the CDF-N , we used the U , B , V , R , / , z! , and HK f pho- 
tometric catalogs from Capak etal. (2004), as well as Spitzer 
IRAC imaging (3.6, 4.5, 5.8, and 8.0^m; Fazio etal. 2004) 
from the GOODS (Dickinson etal., in preparation). For the 
E-CDF-S, we used (1) the 17-bandpass photometry available 
through COMBO- 17, (2) J and K s imaging from the ESO 
Imaging Survey (Olsen etal. 2006), (3) 7, H , and K from 
the MUSYC collaboration (Taylor etal. in preparation; see 
also, Moy etal. 2003; Gawiser etal. 2006), and (4) Spitzer 
IRAC data from the GOODS and SIMPLE teams (Dickin- 
son etal., in preparation; van Dokkum etal., in preparation). 
Using these data, we constructed a rest-frame near-UV-to- 
near-IR SED for each galaxy. For these SEDs, the rest-frame 
U , B, and V filters are well-sampled at all relevant redshifts, 
with the exception of sources at z. ~ 0.7 where the wavelength 
range of the available data is sparse at rest-frame V band. For 
sources at z ~ 0.7, we linearly interpolated our SED to cover 
the V filter. We convolved these SEDs with Johnson U, B, and 
V filter curves and computed rest-frame absolute magnitudes 
for each respective filter following equation 5 of Hogg etal. 
(2002b). For sources in the E-CDF-S, these computed abso- 
lute magnitudes are consistent with those presented by Wolf 
etal. (2004). 

In Figure 1 ( small filled circles ), we present rest-frame 
U-V colors versus My for our galaxies in six redshift ranges. 
For clarity, we also show inset histograms giving the distri- 
bution of rest-frame U - V colors for each redshift interval. 
We utilized the rest-frame U-V color to divide roughly pop- 
ulations of red and blue galaxies. Use of the rest-frame U-V 
color was motivated by Bell et al. (2004a), who note that the U 
and V bandpass pair straddle the 4000 A break, which is par- 
ticularly sensitive to age and metallicity variations of galac- 
tic stellar populations. The dashed lines in Figure 1 show 
the empirically-determined redshift-dependent color divisions 
that separate blue and red galaxy populations; we calculated 
these divisions following § 5 of Bell etal. (2004a): 

(U-V) rest =1.15-0.31z-0.08(M v + 20.7). (1) 

Galaxies having rest-frame U-V less than the values provided 
by equation 1 are often referred to as “blue-cloud” galaxies, 
while those with rest-frame U-V greater than the division are 
called “red- sequence” galaxies. The redshift-dependence of 
the blue-cloud/red-sequence galaxy division is thought to be 
largely due to the evolution of the mean age and dust content 
of the blue-cloud population, with a smaller contribution from 
changes in metallicity. In total, we found 2502 blue-cloud 
galaxies and 757 red-sequence galaxies. 

Generally, the blue-cloud and red-sequence populations are 
composed of late-type and early-type galaxies, respectively 
(see, e.g., Bell etal. 2004b, 2005; McIntosh etal. 2005). 
In support of this point, we have created Figure 2, which 
shows the fraction of E-CDF-S galaxies in our sample hav- 
ing blue-cloud (< dashed histogram ) and red-sequence ( dot- 
ted histogram ) colors as a function of galaxy Sersic index 
n. We utilized the Haussler etal. (2007) Sersic indices, 
which were computed using the GEMS zs 50 images and the 
GALFIT (Peng etal. 2002) two-dimensional light-profile fit- 
ting program. Light-profile studies of large galaxy samples 
have found empirically that a Sersic cutoff of n = 2.5 can 



1 2 3 4 5 6 7 

Sersic Index (n) 


Fig. 2.— Distribution of Sersic indices ( n ) for our initial sample of 
- = 0-1 galaxies in the E-CDF-S having blue-cloud (< dashed histogram) and 
red-sequence ( dotted histogram) colors. The solid vertical line indicates 
n - 2.5, the empirical cutoff between late-type (n £ 2.5) and early-type galax- 
ies (n > 2.5). 


roughly discriminate between late-type and early-type galax- 
ies (e.g., Blanton etal. 2003; Shen etal. 2003; Hogg etal. 
2004). Galaxies with n Zj 2.5 are generally late-type galaxies, 
while the majority of galaxies with n ^ 2.5 are early-types 
( vertical line in Fig. 2). Figure 2 shows that there is reason- 
able agreement between late-type and early-type galaxy pop- 
ulations selected using rest-frame optical colors and Sersic in- 
dices. 

We refined further our division between late-type and early- 
type morphologies by visually inspecting the Zsso-band im- 
ages of our entire sample of 3259 galaxies to see if there were 
obvious cases where the rest-frame optical colors provided an 
inaccurate morphological classification. For example, late- 
type galaxies that are highly inclined to our line of sight may 
experience significant reddening of the young (and blue) disk 
population and will have red-sequence colors. Also, due to 
variations of the stellar populations in galaxies of a given mor- 
phological class, there will be some scatter in rest-frame U-V 
color near the division of the red-sequence and blue-cloud 
regimes. This will lead to a number of “green” galaxies that 
are misclassified morphologically when the classification is 
based solely on rest-frame color. We stress that such a visual 
inspection is only to correct source classifications that were 
obviously misclassified by simply using color-magnitude di- 
agrams, and our resulting sample should not be considered a 
morphologically-selected sample. 

Based on our visual inspection, we found 76 obvious early- 
type galaxies with blue-cloud colors, and 142 obvious late- 
type galaxies with red-sequence colors; we reclassified these 
sources as late-type and early-type galaxies, respectively. We 
inspected histograms of the fraction of blue-cloud and red- 
sequence galaxies that were reclassified as a function of red- 
shift and find no obvious redshift-related trend. In Fig- 
ure 3, we show zsso-band postage-stamp images of ten ob- 
vious early-type galaxies with blue-cloud colors ( top panels ), 
as well as ten obvious late-type galaxies with red-sequence 
colors ( bottom panels). After reclassifying these objects, we 
were left with 2568 late-type galaxies and 691 early-type 
galaxies with z = 0-1 .4. Since we were interested in studying 
the properties of the late- type galaxy populations, we here- 
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Fig. 3.— HST Advanced Camera for Surveys (ACS) zsso example images of ten blue-cloud galaxies that have obvious early-type morphologies ( upper panel ) 
and ten red-sequence galaxies that have obvious late-type morphologies ( lower panel)-, these images demonstrate the aid of the HST coverage. Each postage- 
stamp image has a physical size of 30 kpc x 30 kpc, and in each image we indicate the source name (top), the rest-frame U-V color and redshift of the source 
(lower right), and a vertical bar of length 1"5 for scaling reference (lower left)-, sources harboring an X-ray-detected AGN (see § 4.1) have been annotated with 
the label “AGN” in the lower right. 


after refer to our sample of 2568 late-type galaxies as our main 
sample , which we use in subsequent analyses. 

3. PHYSICAL PROPERTIES AND REDSHIFT EVOLUTION OF 
LATE-TYPE GALAXIES 

The primary goal of this study is to investigate the X-ray 
evolution of normal (i.e., non- AGN) late-type galaxies and 
to determine how this evolution depends upon three intrinsic 
physical properties: optical luminosity, stellar mass, and star- 
formation rate. Below, we describe how we estimated each 
of these physical properties for the galaxies in our main sam- 
ple. We note that the populations of late-type galaxies that 
we are investigating here have been selected via their intrin- 
sic physical properties, which may have changed significantly 
from the observed epoch to the present day. For example, it is 
expected that a significant fraction of the late-type galaxies in 


our main sample will evolve into early-type galaxies via merg- 
ers or passive stellar evolution (e.g., Bell et al. 2007). Ideally, 
we would like to study the evolution of the X-ray properties 
of late-type galaxies while controlling for such changes in the 
physical nature of each galaxy. However, since the details of 
this evolution are highly complex and not well understood for 
a given galaxy, such a task is beyond the scope of this paper. 
We therefore investigate the X-ray evolution of normal late- 
type galaxy populations in relation to their observed intrinsic 
physical properties. 

3.1. Optical Luminosity 

In order to study late-type galaxy samples selected using 
an observable quantity, we made use of the 5-band luminos- 
ity ( L b ). In § 2.3, we computed absolute 5-band magnitudes 
Mb for galaxies in our main sample using photometrically- 




derived SEDs. The B-band emission from a given late-type 
galaxy will be significantly influenced by large populations 
of old ( £ 100 Myr) stars (measured by the stellar mass) as 
well as the younger and less-numerous massive stars that re- 
side in low-obscuration star-forming regions (measured by the 
star-formation rate). In Figure 4 a, we show L B (expressed in 
solar units; where Lb.q = 5.2 x 10 32 ergs s -1 ) versus redshift 
for galaxies in our main sample. The redshift-dependent se- 
lection limit of L b for our sample is set by our zsso < 23 crite- 
rion, and at z = 0.5 (z = 1 .4) this limit corresponds to roughly 
L b « 3 x 10 9 L fl;0 (L fl « 4 x 10 10 For comparison, in 

Figure Aa, we have plotted the values of L B for the Milky 
Way (MW) and local galaxies M101, M82, and the lumi- 
nous infrared galaxy (LIRG) NGC 3256. For the MW, we 
adopted L B = 1.9 x 10 10 L^. 0 as the approximate total disk- 
plus-bulge fi-band luminosity (see Table 15.2 of Gilmore et al. 
1990). For the local galaxies, we utilized optical photometry 
from the Third Reference Catalog of Bright Galaxies (RC3; 
de Vaucouleurs etal. 1991) and adopted distances from the 
IRAS Revised Galaxy Sample (Sanders etal. 2003), which 
gives distances of 6.7, 3.6, and 35.4 Mpc and implies corre- 
sponding B-band luminosities of3.3x 10 10 , 3.8 x 10 9 , and 
4.8 x 10 10 L fi0 for M10 1 , M82, and NGC 3256, respectively. 

3.2. Stellar Mass 

As discussed in § 1 , the X-ray emission from LMXB pop- 
ulations scales with galaxy stellar mass. It is therefore use- 
ful to select late-type galaxies via their stellar masses as a 
means for estimating the LMXB contribution to their X-ray 
emission. To estimate the stellar mass (A/*) of each of our 
galaxies, we exploited the tight correlation between rest- 
frame optical color and stellar mass-to-light ratio (e.g., Bell 
& de Jong 2001; Bell etal. 2003; Kauffmann etal. 2003a; 
see also, Borch etal. 2006). For this calculation, we used a 
combination of rest-frame B-V colors and rest-frame AT-band 
luminosities. X-band luminosities were computed by fitting 
all available optic al/near-IR photometric data to a grid of 276 
synthetic spectra generated by the PEGASE stellar population 
synthesis code (Fioc & Rocca-Volmerange 1997). These tem- 
plates assume a single formation epoch with an exponentially 
decaying star-formation history (time constants r = 1, 2, 4, 
and 7 Gyr) and a Kroupa etal. (1993) IMF; the template grid 
spans ages of 0-15 Gyr. For each galaxy in our main sample, 
we convolved the best-fit template spectrum with the AT-band 
filter function to approximate the rest-frame AT-band luminos- 
ity, L k . This method allowed us to account for significant cur- 
vature in the rest-frame near-IR SEDs of our sources where 
photometric data points did not always overlap with the rest- 
frame X-band. We note that 2546 (^99.1%) of our late-type 
galaxies have at least one photometric data point lying at rest- 
frame wavelengths between 1 and 6 /im, which significantly 
constrains the SED in the rest-frame X-band. We adopted the 
prescription outlined in Appendix 2 of Bell etal. (2003) to 
estimate Af* using rest-frame B-V color and X-band lumi- 
nosity: 

\ogM+/M Q = \ogL K /L K . Q +0A35(B-V)-0.306. (2) 

The numerical constants in equation 2 were supplied by Ta- 
ble 7 of Bell etal. (2003) and are appropriate for our choice 
of B-V color and L K \ the normalization has been adjusted 
by -0.1 dex to account for our adopted Kroupa (2001) IMF 
(see § 1). We compared our stellar-mass estimates with those 
computed by Borch et al. (2006) for 1758 (^68% of our main 
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Fig. 4.— Rest-frame B-band luminosity Lb ( a ), stellar mass A/* ( b ), and 
star-formation rate SFR (c) versus redshift for sources in our main sample. 
For comparison, in each panel we have indicated the properties of the MW 
and local galaxies M101 , M82, and NGC 3256. In Figure 4c, we have in- 
cluded only SFRs for the 888 late-type galaxies in our main sample that have 
24/im counterparts (see § 3.3 for details); the dashed curve indicates roughly 
the SFR detection limit for sources without 24/im counterparts. X-ray- 
detected normal galaxies and AGNs have been indicated with open circles 
and diamonds, respectively. The thick gray rectangles indicate regions where 
galaxy populations were selected for X-ray stacking (see §§4.2 and 5.1). For 
each galaxy. M* and SFR were computed following equations 2 and 3, re- 
spectively. 


sample) galaxies in the COMBO- 17 survey and find excel- 
lent agreement between methods, with an overall scatter of 
^0.2 dex. 

In Figure 4 b, we show Af* versus redshift for our sample of 
late- type galaxies. We note that our stellar-mass estimates are 
broadly limited by the zgso < 23 criterion used in our sample 
selection. At z- 0.5 (z = 1 .4) our samples are representative 
for late-type galaxies with Af* £ 10 9 Af 0 (Af* £ 10 10 Af o ). For 
comparison, in Figure 4b we have plotted stellar masses for 
the MW and local galaxies M101 , M82, and NGC 3256. We 
adopted a stellar mass of 5 x 10 10 Af 0 for the MW (e.g., Ham- 
mer et al. 2007). For the local galaxies, we estimated Af* using 
a similar relation to equation 2, except using B-V color and 
fl-band luminosity (see Appendix 2 of Bell etal. 2003 for de- 
tails); for this computation, we utilized B-V colors from RC3 
and values of L B determined in § 3.1 . This computation gives 
stellar masses of 1 .3 x 10 10 , 5.4 x 10 9 , and 7.8 x 10 10 Af 0 for 
M101 , M82, and NGC 3256, respectively. 

3.3. Star- Formation Rates 

Since the X-ray emission from normal late- type galaxies 
is known to be strongly correlated with SFR, it is of partic- 
ular interest to understand how changes in SFR have con- 
tributed to the X-ray evolution of the normal late-type galaxy 
population. To calculate SFRs for the galaxies in our main 
sample, we utilized estimates of both the dust-uncorrected ul- 
traviolet luminosities (L uv ) originating from young stars and 
the infrared luminosities (8-1000/im; Lir) from the dust that 



obscures UV light in star-forming regions (see, e.g., Ken- 
nicutt 1998 for a review). The former quantity was com- 
puted following Luv = 3.3^^(2800 A), where l „{ 2800 A) is 
the rest-frame 2800 A monochromatic luminosity (see § 3.2 
of Bell etal. 2005). ^(2800 A) was approximated using 
the template SEDs that were constructed in § 3.2. The lat- 
ter quantity (Lir) was computed using observed- frame 24^m 
flux densities (i.e., rest-frame 24/im/l +z) from observations 
with the MIPS (Rieke etal. 2004) camera onboard Spitzer. 
We note that generally the dust-obscured star-formation ac- 
tivity, probed here using Lir, can be measured using either 
the UV spectral slope or optical nebular recombination lines 
(e.g., Balmer emission lines and O II); however, the data avail- 
able for our z ~ 0-1 .4 late-type galaxies lack the rest-frame 
UV and spectroscopic coverage needed to provide such useful 
measurements. 

Over the CDFs, deep 24pm observations were available for 
the GOODS fields lim ~ 30 ply, 6a). 13 For the remain- 

ing area (covering the outer regions of the E-CDF-S), shal- 
lower 24pm observations (/24^ m , lim ~ 120 ply) were avail- 
able through SWIRE (Lonsdale etal. 2003). We matched the 
positions of the 2568 late-type galaxies in our main sample 
with those from the available 24pm source catalogs, requir- 
ing that the HST and Spitzer centroids be offset by no more 
than lf'5. We found successful matches for 888 (^35%) of 
our late-type galaxies. For these sources, we converted 24^m 
flux densities to Lir following the methods outlined in Pa- 
povich & Bell (2002). Briefly, we utilized the entire grid of 
64 infrared SEDs provided by Dale etal. (2001) to estimate 
the mean conversion factor rj(z ), which transforms observed- 
frame 24jLtm luminosity ul u (24pm/\+z) to Lir as a function 
of redshift [i.e., Lir = p(j)vl u (24 pm j\ + z)]. The mean con- 
version factor spans the tight range of 77 = 7.2-12.2 over the 
entire redshift range z = 0-1 .4. We note that different choices 
of infrared SEDs yield similar results at z £ 1 , but become 
significantly discrepant at z £ 1 (see, e.g., Fig. 2 of Papovich 
& Bell 2002). 

Using our estimates of Luv and Lir, we calculated star- 
formation rates for galaxies in our sample using the following 
equation: 

SFR(M 0 yr') = 9.8 x l<T n aiR+Luv), (3) 

where Lir and Luv are expressed in units of the solar bolo- 
metric luminosity (L 0 = 3.9 x 10 33 ergs s -1 ). Equation 3 
was adopted from § 3.2 of Bell etal. (2005) and was de- 
rived using PEGASE stellar-population models, which as- 
sumed a 100 Myr old population with constant SFR and a 
Kroupa (2001) IMF (see Bell 2003 for further details). We 
note that for the majority of our galaxies Lir > Luv , and for 
sources with z < 0.5 and z > 0.5, the median IR-to-UV ratio 
Ur/L\jv = 2.8 and 7.1, respectively. Our z < 0.5 galaxies have 
median Lir = 1.1 x 1O 1O L 0 , which is characteristic of objects 
that are found in the nearby D < 10 Mpc universe. By con- 
trast, our z > 0.5 galaxies have median Lir = 9.3 x 10 
similar to the local LIRG population found at D > 10 Mpc. 
In Figure 4c, we show the distribution of SFRs for the 888 
galaxies in our main sample that had 24^m counterparts. We 
note that for sources in our main sample that were within 
the GOODS regions, where the 24pm observations are most 
sensitive, the infrared detection fraction drops from ^100% 
for galaxies with zsso = 20 ±0.2 to «20% for galaxies with 

13 Available at http://ssc.spitzer.caltech.edu/legacy/goodshistory.html. 


Z850 = 22.8 ±0.2. This demonstrates that our SFR complete- 
ness is limited primarily by the 24pm sensitivity limit and 
that our sample of zgso < 23 late-type galaxies is highly rep- 
resentative of galaxy populations above the apparent redshift- 
dependent SFR limit shown in Figure 4c (< dashed curve). At 
z = 0.5 (z = 1 .4), this limit corresponds to SFR ^ 1 Mq yr -1 
(SFR ^ 15 Mq yr -1 ). For comparison, in Figure 4c we 
have plotted the SFRs for the MW and local galaxies M101, 
M82, and NGC 3256. The SFR of the MW was taken to be 
«2 Mq yr -1 (McKee & Williams 1997), and the SFRs of the 
local galaxies were computed using equation 3. Values of 
2800 A) were approximated from SEDs available through 
the NASA/IPAC Extragalactic Database (NED) 14 and L JR was 
taken from the IRAS Revised Bright Galaxy Sample (Sanders 
etal. 2003). We find SFR = 2.5, 6.2, and 32.2 M 0 yr 1 for 
M101,M82, and NGC 3256, respectively. 

As a consistency check on our UV-plus-infrared SFR es- 
timates SFR(UV+IR), we calculated radio-derived SFRs us- 
ing 1.4 GHz observations, SFR(1.4 GHz), following equa- 
tion 7 of Schmitt etal. (2006), which we adjusted to be con- 
sistent with our adopted Kroupa (2001) IMF. We matched 
sources in our main sample to 1.4 GHz catalogs, which 
were derived from observations using the VLA in the CDF-N 
(«30 ply, Richards etal. 1998) and the ATCA in the 
E-CDF-S (^60 pJy; Afonso etal. 2006; Rovilos etal. 2007). 
Using a matching radius of 1"5, we found a total of 54 
radio sources coincident with our late-type galaxies. We 
found that 48 (?^89%) of the radio-detected sources had 
24pm counterparts, allowing for reasonable comparison be- 
tween derived SFRs. For these sources we found rea- 
sonable agreement between SFRs derived from UV-plus- 
infrared and radio measurements, with a mean ratio of 
SFR( 1 .4 GHz)/SFR(U V +IR) = 1.5 ± 1.0. A large number 
of late-type galaxies (840 sources) in our main sample were 
detected in the 24pm observations that were not detected at 
1.4 GHz. All but three of these sources had SFR(1.4 GHz) 
upper limits that were consistent with that expected from es- 
timates of SFR(UV+IR) (see § 4.1 for further details). 

4. ANALYSIS 

4.1 . X-ray-Detected Late-Type Galaxies and AGN 
Identification 

We utilized the multiwavelength observations in the CDFs 
to obtain a census of the active galaxies in our main sample. 
We began by matching the optical positions of our galaxies to 
the X-ray positions of point sources in the CDF catalogs of 
Alexander etal. (2003) 15 for the ^2 Ms CDF-N and «1 Ms 
CDF-S and Lehmer et al. (2005b) 16 for the ^250 ks E-CDF-S. 
For a successful match, we required that the optical and X-ray 
centroids be displaced by no more than 1 .5 times the radius of 
the Chandra positional error circles (80%-90% confidence), 
which are provided in each respective catalog. We note that 
for a small number of the galaxies in our sample at z 0.3, 
moderately luminous off-nuclear X-ray sources (e.g., ultralu- 
minous X-ray sources [ULXs]) that are intrinsically related to 
the galaxies may lie outside of our adopted matching radius; 
we utilized the off-nuclear X-ray source catalog of Lehmer 

14 Available at http://nedvvvvw.ipac.caltech.edu/. 

15 See http://wwvv.astro.psu.edu/users/niel/hdf/hdf-chandra.html for the 
relevant source catalogs and data products for the zz2 Ms CDF-N and ssl Ms 
CDF-S. 

16 See http://www.astro.psu.edu/users/niel/ecdfs/ecdfs-chandra.html for 
the relevant source catalogs and data products for the ^250 ks E-CDF-S. 
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FIG. 5.— Count-rate ratio in the 2-8 keV to 0.5-2 keV band- 
passes ($ 2-8 kcv/^ 0 . 5-2 kev) versus the logarithm of the 0.5-8 keV flux 
(log/0.5-8 kev) for X-ray-4etected sources in our main sample. Off-nuclear 
X-ray sources catalogued by Lehmer et al. (2006) have been highlighted with 
open squares. The shaded region represents sources with effective photon 
indices T e ff ± 1 ; we classified these sources as AGN candidates (see § 4, cri- 
terion 1). For reference, we have plotted lines corresponding to r e ff = 0.5, 1 , 
and 2 {dashed lines) . 



etal. (2006) to identify such galaxies and assign X-ray prop- 
erties (see details below). 

The Chandra source catalogs were generated using 
wavdetect (Freeman etal. 2002) with false-positive prob- 
ability thresholds of 1 x 10 -7 and 1 x 10 -6 for the Alexan- 
der etal. (2003) and Lehmer etal. (2005b) point-source cat- 
alogs, respectively. However, as demonstrated in § 3.4.2 of 
Alexander etal. (2003) and § 3.3.2 of Lehmer etal. (2005b), 
legitimate lower significance X-ray sources, detected by run- 
ning wavdetect at a false-positive probability threshold of 
1 x 10~ 5 , can be isolated by matching with relatively bright 
optical sources; therefore, when matching our late-type galax- 
ies to X-ray-detected sources, we utilized this technique. The 
sky surface density for late-type galaxies in our main sample 
ranges from «23000 deg~ 2 in the CDF-N to ^8800 deg"Mn 
the «250 ks E-CDF-S. The large difference between these 
source densities is primarily due to differences in applied red- 
shift cuts (see § 2.2). Given the fact that the positional un- 
certainties are generally small ( £ 1 "5) for sources within 
6f0 of the Chandra aim points, as is the case for sources 
in our main sample, the corresponding estimated number of 
spurious matches is small. We estimate that when using 
wavdetect with a false-positive probability threshold of 
1 x 10" 5 to search for sources in three Chandra bandpasses 
(0.5-2 keV, 2-8 keV, and 0.5-8 keV), we expect ?^1 .8 spuri- 
ous matches. When including the off-nuclear sources from 
Lehmer etal. (2006), we expect an additional ^0.5 false 
sources; this brings our total spurious matching estimate to 
~2.3 sources for our main sample. 

Using the matching criteria above, we find that 225 late- 
type galaxies are detected in at least one of the 0.5-2 keV, 
2-8 keV, or 0.5-8 keV bandpasses. Out of these 225 galaxies, 
12 are known off-nuclear X-ray sources from Lehmer etal. 
(2006). Since only one off-nuclear source is detected for each 
host galaxy, we assume that the off-nuclear point-source dom- 
inates the total X-ray emission from each host galaxy. We 
therefore adopted the X-ray properties presented in Table 1 of 
Lehmer et al. (2006) for each off-nuclear host galaxy. 

The unprecedented depths of the CDFs allow for the indi- 


vidual X-ray detection of luminous normal galaxies over the 
entire redshift range of our main sample (z = 0-1 .4); however, 
the majority of the X-ray-detected sources in even the most- 
sensitive regions of the CDFs will be distant (z ^ 0.5) AGNs 
(e.g., Bauer etal. 2004), which we want to separate from our 
normal late-type galaxy sample. We identified AGN candi- 
dates using four primary criteria, which utilize (1) X-ray hard- 
ness to identify luminous obscured sources, (2) X-ray-to- 
optical flux ratios to identify additional relatively-unobscured 
AGNs, (3) the X-ray-SFR correlation to identify additional 
lower-luminosity AGNs that are significantly influencing the 
total X-ray emission, and (4) the combination of radio and in- 
frared properties to identify additional Compton-thick AGNs 
and radio-loud AGNs. As a final check on our AGN identifi- 
cations we utilized optical spectroscopic information to iden- 
tify sources with obvious AGN signatures (see criterion 5 be- 
low). In the sections below we provide details of each criteria. 

1. X-ray Hardness.— One unique signature of moderately 
obscured (Ah ^ 10 22 cm -2 ) AGN activity is a hard X-ray 
spectrum. For normal galaxies, the collective emission from 
X-ray binaries dominates the total 0.5-8 keV power output. 
On average, these sources have observed power-law X-ray 
SEDs with spectral slopes of T « 1 .5-1 .7 for LMXBs (e.g., 
Church & Balucinska-Church 2001; Irwin etal. 2003) and 
T « 1-2 for HMXBs and ULXs (e.g., Sasaki etal. 2003; Liu 
& Mirabel 2005; Liu etal. 2006); the presence of a signifi- 
cant hot interstellar gas component will steepen the resulting 
X-ray spectral slope (i.e., produce larger r e ff). To identify ob- 
scured AGNs in our sample effectively, we flagged sources 
having effective photon indices of r c ff ^ 1 as AGN candi- 
dates. We determined T e ff using the 2-8 keV to 0.5-2 keV 
hardness ratio $ 2 -8 kev/ $ 0 . 5-2 kev, where $ is the count rate 
for each bandpass. A few sources have only 0.5-8 keV detec- 
tions. Since these sources were not detected in the 0.5-2 keV 
bandpass, our most sensitive bandpass, there must be a sig- 
nificant contribution from the 2-8 keV bandpass such that 
$2-8 kev/$0 .5—2 keV 1 (Faff <0.9). We therefore classified 
these sources as AGN candidates. In Figure 5, we show the 
$ 2-8 kev/$o. 5-2 kev hardness ratio versus the logarithm of the 
0.5-8 keV flux log/o. 5 - 8 kev for the X-ray-detected sources 
in our main sample. The shaded area highlights the region 
corresponding to r e ff £ 1 . Sources that had only upper lim- 
its on $ 2-8 kev / $ 0 . 5-2 kev that lie in the shaded region were 
not classified as AGN candidates, while those with lower lim- 
its on $ 2-8 kev / $ 0 . 5-2 kev that were in the shaded region were 
flagged as likely AGNs. We note that occasionally an X-ray 
luminous ULX that is too close to the nucleus of its host 
galaxy to be identified as off-nuclear may have r e ff < 1 and 
therefore be classified as an AGN candidate via this criterion; 
however, since only one identified off-nuclear X-ray source 
has r e ff < 1 ( open squares in Fig. 5; n.b., upper limits are not 
included), we do not expect such rare sources to significantly 
affect our results. Using criterion 1 , we identified a total of 7 1 
obscured AGN candidates. 

2. X-ray-to-Optical Flux Ratio.— Detailed analyses of the 
X-ray spectra of luminous AGNs in the Ms CDF-S show 
that the intrinsic AGN power-law photon index is relatively 
steep, (T) = 1.75 ±0.02 (e.g., Tozzi etal. 2006). Therefore, 
luminous AGNs having column densities of Ah £ 10 22 cm -2 
will often have effective photon indices of r c ff > 1 and would 
not have been classified as potential AGNs by criterion 1. 
In order to identify luminous AGNs with Ah ^ 10 22 cm -2 , 
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FIG. 6.— /?-band magnitude versus log/ 0 . 5-8 keV for X-ray-detected 
sources in our main sample. Open circles represent sources that were classi- 
fied as AGN candidates by criterion 1 (see § 4.1 and Fig. 5), and filled circles 
represent all other sources; open squares have the same meaning as in Fig- 
ure 5. The shaded area represents the region where log(/o. 5_8 kev///?) > -1; 
we classified sources in this region as AGN candidates (see §4.1, criterion 2). 
For reference, the dashed lines represent log(/o. 5-8 kev//?) = -2, -1 , and 1 . 

we utilized the X-ray-to-optical flux ratio (/ 0 . 5-8 kev / Ir) as 
a discriminator of AGN activity (e.g., Maccacaro etal. 1988; 
Homschemeier etal. 2000; Bauer etal. 2004). We identified 
sources with log(/o. 5-8 kev ///?)> -1 (see criterion 3 below for 
justification) as unobscured AGN candidates. In Figure 6, we 
show the fl-band magnitude versus log/ 0 . 5-8 kev for sources 
in our main sample; the shaded area shows the region where 
log(/o.5-8 k cv//z?) > -1 • Sources that were classified as AGN 
candidates via criterion 1 are denoted with open symbols. In 
total, 43 X-ray-detected sources satisfied criterion 2, and 16 
of these sources were uniquely identified using this criterion 
(i.e., not identified by criterion 1). 

3. X-ray-to-SFR Correlation.— Taken together, criteria 1 
and 2 provide an effective means for identifying AGNs that 
are affected by large absorption column densities (criterion 1) 
and those that are notably X-ray overluminous for a given 
optical luminosity (criterion 2). However, these criteria will 
still miss moderately luminous unobscured AGNs that have 
log(/o. 5-8 k ev//z?) < -1 (see, e.g., Peterson etal. 2006). Al- 
though an accurate classification for all such sources is cur- 
rently not possible, the situation can be mitigated using the 
available multi wavelength data. We therefore exploited the 
correlation between X-ray luminosity Lx and SFR (see § 1) to 
identify additional AGN candidates in our main sample that 
have significant X-ray excesses over what is expected based 
on the L x -SFR correlation. In order to calculate the rest-frame 
luminosity ( Le k -e 2 \ where £j and £2 are the photon-energy 
lower and upper bounds, respectively) of a source having a 
power-law SED, we used the following equation: 

Le { -e 2 - ^d±fE\-E 2 (\ +/) r-2 > ( 4 ) 

where fE { -E 2 is the observed- frame emission in the E 1 -E 2 
bandpass and cIl is the luminosity distance. Using equa- 
tion 4 and an adopted photon index of T = 2, we calcu- 
lated 0.5-8 keV luminosities for X-ray-detected sources in 
our sample. In Figure 7, we show the logarithm of the 
0.5-8 keV luminosity logLo. 5-8 kev versus SFR (computed 
following equation 3) for galaxies in our main sample that 
had 24^m counterparts (see § 3.3). Several estimates of the 
Lx -SFR correlation have been shown for reference (Bauer 
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Fig. 7.— Logarithm of the 0.5-8 keV luminosity Lo. 5-8 kev versus SFR for 
X-ray-detected sources in our sample that have 24//m counterparts. Open 
circles represent sources that were classified as AGN candidates via crite- 
ria 1 and 2 (see § 4.1 and Figs. 5 and 6), and filled circles represent all 
other sources; open squares have the same meaning as in Figure 5. We have 
shown the X-ray-SFR relations calibrated by Bauer et al. (2002) , Ranalli et al. 
(2003), Gilfanov etal. (2004a), and Persic & Rephaeli (2007); each respec- 
tive curve has been annotated in the figure. The shaded area above the thick 
dashed line represents the region where Lo. 5-8 kev is three times larger than 
its value predicted by Persic & Rephaeli (2007); we classified 0.5-8 keV de- 
tections that lie in this region as AGN candidates (see § 4.1 , criterion 3). 


etal. 2002a; Ranalli etal. 2003; Gilfanov etal. 2004a; Per- 
sic & Rephaeli 2007 , hereafter PR07); these correlations have 
been corrected for differences in X-ray bandpass and SED as 
well as adopted IMFs. Hereafter, we adopt the Lx-SFR corre- 
lation from PR07 for comparisons; however, the use of other 
Lx-SFR correlations would yield similar results and conclu- 
sions. Open squares show the locations of the galaxies hosting 
off-nuclear sources, which appear to be preferentially located 
near the L X -SFR correlation. Open circles indicate sources 
that were identified as AGN candidates via criteria 1 and 2. 
Generally, these AGNs have Lo.5-8 kcv/SFR £ 3 times that 
predicted by the PR07 L X -SFR correlation ( thick dashed line), 
which corresponds to a factor of >,2.5 times the RMS scat- 
ter of the PR07 Lx-SFR correlation. We therefore classified 
all sources in this regime {shaded region) as AGN candidates. 
We note that sources having only 0.5-8 keV upper limits were 
not classified as AGNs. However, sources that were detected 
in the 0.5-8 keV bandpass that had only upper limits on the 
SFR were classified as AGN candidates if they were within 
the shaded region of Figure 7. Using this criterion, we identi- 
fied 101 potential AGN candidates, of which 33 were unique 
to criterion 3 . For the 1 24 X-ray-detected sources that were 
not classified as AGNs, we found that 19 sources («15%) had 
X-ray spectral properties indicative of AGNs (criterion 1) and 
only one source (^0.8%) had an X-ray-to-optical flux ratio 
elevated in comparison to normal galaxies (criterion 2). 

4. Infrared and Radio Properties .— In addition to the 
X-ray-detected sources, we also expect there to be additional 
AGNs that lie below the X-ray detection threshold that are not 
identified here. In § 3.3, we noted using 48 sources that the 
SFRs derived from the UV-plus-infrared emission are gener- 
ally consistent with those derived from radio (1 .4 GHz) obser- 
vations; however, there are a few sources that are significantly 
scattered (i.e., by more than a factor of five times the intrin- 
sic scatter) outside of the SFR(UV+IR)-SFR(1 .4 GHz) cor- 
relation. These sources have either (1) excess infrared emis- 



Fig. 8.— Logarithm of the 0.5-8 keV luminosity Lo 5 _g kev versus redshift 
for the 225 X-ray-detected sources in our main sample. Open circles indicate 
all of the AGN candidates that were isolated using the four criteria outlined 
in § 4.1, and filled circles represent potential normal galaxies; open squares 
have the same meaning as in Figure 5. The estimated 0.5-8 keV detection 
limits for the «2 Ms CDF-N, the « 1 Ms CDF-S, and the &250 ks E-CDF-S 
have been indicated with dotted, dashed, and dot-dashed curves, respectively. 


sion due to reprocessed dust emission from a highly-obscured 
(Nu 10 24 cm -2 ) moderately luminous (intrinsic Lx ~ 10 42 - 
10 43 ergs s~ l ) Compton-thick AGN (e.g.,Daddi etal. 2007 a, b) 
or (2) an excess of radio emission due to prominent radio- 
emitting jets (e.g., Snellen & Best 2001). The reasonable 
agreement between SFR(UV+IR) and SFR(1 .4 GHz) suggests 
that such AGNs are not prevalent in our late-type galaxy sam- 
ple. Ideally we would also like to explore how UV-corrected 
SFRs compare to the UV-plus-infrared SFR to search for po- 
tential Compton-thick AGNs following the approach used by 
Daddi etal. (2007a, b) for BzK galaxies at 2; however, 
since our galaxies generally lie at much lower redshifts that 
the Daddi etal. (2007a, b) sources, we are unable to constrain 
well UV spectral slopes with our available photometry, which 
would be required to make dust-corrections to the observed 
UV SFRs (see § 3.3 for discussion). For our sample of 24/rni- 
detected galaxies, we identified three radio-excess galaxies 
with SFR(1 .4 GHz)/SFR(UV+IR) > 5 and three IR-luminous 
galaxies with SFR(1.4 GHz)/SFR(UV+IR) <1/5. Using 
these criteria we identified three new AGN candidates, two ra- 
dio luminous (J33 146.6-275735 and J33244.3-275141) and 
one IR luminous (J3 3240 .0-2742 14); one of these sources, 
the radio luminous source J33244.3-275141, was detected in 
the 0.5-2 keV and 0.5-8 keV bandpasses. 

As an additional test, we utilized IRAC photometry (see 
§ 2.3) to search for infrared power-law sources having near- 
IR spectral properties characteristic of luminous AGNs (e.g., 
Alonso-Herrero etal. 2006; Donley etal. 2007). When 
searching for power-law sources, we adopted the criteria dis- 
cussed in Donley etal. (2007). We found that no sources in 
our sample satisfied these criteria, which is consistent with the 
finding that most IRAC power-law sources reside at z ~ 1 . 

5. Optical Spectroscopy.— As a final check on our AGN 
classifications, we searched the optical spectroscopic catalogs 
available for CDF sources (e.g., Barger etal. 2003; Le Fevre 
etal. 2004; Szokoly etal. 2004; Wirth etal. 2004; Vanzella 
etal. 2005,2006) to isolate additional luminous AGNs in our 
sample. In total, we found 15 galaxies in our main sam- 
ple that were classified as AGNs via optical spectroscopy 


and all of these sources had been identified as AGN candi- 
dates by the previous criteria. We note that the majority of 
the X-ray-detected AGNs have moderate luminosities (intrin- 
sic Lx ~ 10 41 — 10 43 ergs s -1 ) and therefore often have high- 
excitation AGN emission lines that are too faint with respect 
to stellar emission to be identified via optical spectroscopy 
(e.g., Moran etal. 2002). Therefore, it is not surprising that 
we do not find any additional AGNs using this criterion. 

To summarize , in the X-ray band we have detected a total 
of 225 («9%) late-type galaxies out of the 2568 sources in 
our main sample. Using the criteria presented above, we clas- 
sified 121 X-ray-detected sources as AGN candidates (with 
an additional two X-ray-undetected AGN candidates via cri- 
terion 4). The remaining 104 X-ray-detected sources that we 
do not classify as AGN candidates are considered to be nor- 
mal late-type galaxies, and we include these galaxies in sub- 
sequent X-ray stacking analyses (see details in § 5 below). 
Thus we use 2447 late-type galaxies in our stacking analyses. 
In Figure 8, we show logLo. 5 _ 8 keV versus redshift for X-ray- 
detected sources in our main sample, and in Table 1 {available 
electronically) we summarize their properties. AGN candi- 
dates are denoted with open circles and normal galaxies are 
plotted with filled circles. We note that the above criteria 
are not completely sufficient to classify all X-ray-detected 
sources that are truly AGNs as AGN candidates (see, e.g., 
Peterson etal. 2006). Such a misclassification is possible 
for low-luminosity AGNs that are only detected in the more 
sensitive 0.5-2 keV bandpass, which have 2-8 keV emission 
too weak for an accurate classification. In § 5.2, we use the 
2-8 ke V AGN fraction as a function of X-ray luminosity to ar- 
gue quantitatively that we do not expect misclassified AGNs 
(detected only in the 0.5-2 keV bandpass) and low-luminosity 
AGNs below the X-ray detection limit to have a serious im- 
pact on our results. 

4.2. X-ray Stacking Analyses of Normal Late -Type Galaxy 
Populations 

The majority of the normal late-type galaxies that make up 
our main sample were not detected individually in the X-ray 
bandpass. In order to study the mean X-ray properties of 
these sources, we implemented stacking analyses of galaxy 
populations selected by their physical properties (see § 3). 
We divided our main sample into subsamples (to be used for 
stacking) of normal late-type galaxies selected by both physi- 
cal properties (i.e., 5-band luminosity, stellar mass, and star- 
formation rate) and redshifts. In Figure 4, we have highlighted 
the divisions of our sample with thick gray rectangles, and for 
normal late-type galaxies in each region, we used X-ray stack- 
ing analyses to constrain average properties. 

For each of the subsamples defined above, we performed 
X-ray stacking in each of the three standard bandpasses (i.e., 
SB, HB, and FB; see § 1). We expect these bandpasses 
to sample effectively power-law X-ray emission originating 
from X-ray binaries (i.e., HMXBs and LMXBs) with a minor 
contribution from hot interstellar gas in the SB for late-type 
galaxies at z ~ 0.5. In our analyses, we used data products 
presented in Alexander etal. (2003) for the ^2 Ms CDF-N 
and « 1 Ms CDF-S and Lehmer et al. (2005b) for the ^250 ks 
E-CDF-S (see footnotes 15 and 16). Our stacking procedure 
itself was similar to that discussed in § 3 of Steffen etal. 
(2007). This procedure differs from past stacking analyses 
(e.g., Lehmer etal. 2007) in how the local X-ray background 
of each stacked sample is determined, and produces results 







Table 1 . X-ray Detected Late-Type Galaxies: Source Properties 


Source Name zsso Lb M * SFR £ 0 . 5-8 kcv / 0 . 5-8 kev 

(J2000.0) (mag) z (log Lb q) (logM©) (M© yr -1 ) Survey (ks) (log ergs cm -2 s _1 ) 

(1) (2) (3) (4) (5) (6) (7) (8) (9) 


J033 1 32 .8 1 -280 1 1 5 .9 17.19 0.15P 10.39 10.49 <5.89 E-CDF-S 03 217 -15.06 

J033 132.84-280107 .5 22.17 0.15P 8.85 8.44 0.23 E-CDF-S 03 217 -15.45 

J033 139. 20-280222. 3 17.06 0.25P 10.94 11.07 1.40 E-CDF-S 03 206 <-15.15 

J033 139 .97-274 157.0 20.14 0.47 s 10.14 10.43 <19.93 E-CDF-S 02 226 -13.46 

J033 14 1.79-275635 .2 19.58 0.5 IP 10.84 10.54 <13.17 E-CDF-S 03 239 -15.02 


NOTE. — Col.(l): Chandra source name. Col.(2): ACS zsso-band magnitude. Col.(3): Rcdshift estimate. Superscripts “s” and “p” indicate spectroscopic and photometric redshifts, 
respectively (see § 2 for details). Col.(4): Logarithm of the rest-frame 5-band luminosity in units of Lb.q . Col.(5): Logarithm of the stellar mass in units of Af© . Col.(6): Star-formation 
rate in units of Mq yr -1 . Col .(7): Survey field in which each source was identified. For E-CDF-S identifications, the associated field number (i.e., 01-04) indicates the Chandra pointing 
within which the source was detected (see Lehmer et al. 2005b for details). Col.(8): Effective 0.5-8 keV exposure time (in units of ks). Col.(9)-(l 1): Logarithm of the 0.5-8 keV, 
0.5-2 keV, and 2-8 keV flux in units of ergs cm -2 s -1 . Col .( 1 2)— ( 14): Logarithm of the 0.5-8 keV, 0.5-2 keV, and 2-8 keV rest-frame luminosity in units of ergs s -1 . Col.(15): 
Effective photon index (r c f f ). Col.(16): Logarithm of the 0.5-8.0 keV to 5-band flux ratio. Col.(17): AGN candidate (Y/N)? Col.(18): AGN selection criteria used (i.e., 1-5; see § 4.1 
for details). All 225 entries of Table 1 are available in the electronic edition. A portion is shown here for guidance regarding its form and content. 


Table 2. Stacked Late-Type Normal Galaxies: Mean Properties 


Selection Type 
0) 

£mean 

(2) 

Ngal 

(3) 

Ndet 

(4) 

Net Source Counts (S-B) 

0.5-8 keV 0.5-2 keV 2-8 keV 
(5) (6) (7) 

Signal-to-Noise (S/N) 

0.5-8 keV 0.5-2 keV 2-8 keV 
(8) (9) (10) 

log Lb/Lb.q =9.5-10.0 

0.23 ±0.06 

120 

9 

213.3 

158.7 

54.9 

8.6 

9.6 

3.0 

log L b /Lb.q =9.5-10.0 

0.46 ±0.06 

363 

1 

176.7 

152.6 

25.3 

4.8 

6.8 

0.9 

log Lb! Lb.q = 10.0-10.5 

0.24 ±0.02 

28 

2 

66.9 

46.5 

20.2 

5.2 

5.3 

2.1 

log Lb/Lb,q = 10.0-10.5 

0.33 ±0.03 

33 

3 

82.9 

68.5 

14.9 

5.9 

6.7 

1.5 

log Lb/Lb, © = 10.0-10.5 

0.41 ±0.03 

67 

7 

204.2 

144.7 

59.8 

9.3 

9.7 

3.8 

logL B /L B . e = 10.0-10.5 

0.50 ±0.02 

158 

6 

347.5 

240.0 

109.1 

10.7 

11.3 

4.4 

log/.g/Z.B, 0 = 10.0-10.5 

0.57 ±0.02 

100 

3 

128.0 

109.0 

18.7 

6.0 

7.5 

1.2 

\ogL B /L B .@ = 10.0-10.5 

0.67 ±0.02 

71 

2 

154.1 

93.5 

60.6 

6.3 

6.1 

3.2 

\o%L b /L b ,q = 10.0-10.5 

0.75 ±0.02 

62 

2 

156.3 

107.4 

48.8 

6.4 

6.9 

2.6 

\o%L b /L b ,q = 10.0-10.5 

0.84 ±0.02 

57 

3 

105.9 

64.8 

41.2 

4.4 

4.5 

2.2 


Note. — Col.(l): Physical property (i.e., Lb, M* , and SFR) used to select the stacked sample. Col.(2): Mean redshift (Zmean) and 1 a standard deviation. Col.(3): Number of 
galaxies stacked (N g . d[ ). Col.(4): Number of X-ray-detected normal galaxies stacked. Col.(5)-(7): Net source counts ( S-B ) for the 0.5-8 keV, 0.5-2 keV, and 2-8 keV bandpasses. 
Col.(8H10): Signal-to-noise ratio (S/N) for the 0.5-8 keV, 0.5-2 keV, and 2-8 keV bandpasses. Col.(l 1 )— (1 3): Logarithm of the mean 0.5-8 kcV, 0.5-2 kcV, and 2-8 keV flux 
in units of ergs cm -2 s _l . Col.(14)-(16): Fraction of the mean 0.5-8 keV, 0.5-2 keV, and 2-8 keV flux originating from the X-ray-undetected galaxies. Col.(17)— (19): Logarithm 
of the mean 0.5-8 keV, 0.5-2 keV, and 2-8 keV rest-frame luminosity. Col.(20): Mean effective photon index (Tefr). Col.(21): Logarithm of the 0.5-8.0 keV to 5-band flux ratio. 
Col.(22)-(24): Mean values of L B , M* , and SFR for samples selected by L fl ,M* , and SFR, respectively. Col.(25)-(27): Logarithm of the ratios L\ JL B (ergs s' 1 L~ B q),L\/M * (ergs s -1 
M~q), and L x /SFR [eigs s -1 (Af© yr -1 )' 1 ] for samples selected by L a , M* , and SFR, respectively. Col.(28)-(30): Estimated fraction of the mean 0.5-8 keV, 0.5-2 keV, and 2-8 keV 
stacked emission originating from undetected AGNs. All 44 entries of Table 2 are available in the electronic edition. A portion is shown here for guidance regarding its form and 
content. 


that are in good agreement with the method discussed in § 2.2 
of Lehmer etal. (2007). For completeness, we have outlined 
this procedure below. 

Using a circular aperture with radius R ap = 1"5, we ex- 
tracted Chandra source-plus-background counts 5, and expo- 
sure times 7} (in units of cm 2 s) for each galaxy using images 
and exposure maps, respectively. For a given source, we used 
only Chandra pointings with aimpoints (see footnote 1 2) that 
were offset from the source position by less than 6'0; here- 
after, we refer to this maximum offset as the inclusion radius, 
T?inci- 17 Since the Chandra PSF increases in size with off- 
axis angle and degrades the sensitivity for sources that are 
far off-axis, our choices of R ap and R[ nc \ have been chosen 
to give the maximal stacked signal with the majority of the 
PSFs being sampled by our stacking aperture (see § 2.2 and 
Fig. 3 of Lehmer etal. 2007 for justification). We estimate, 
based on stacked images of z ~ 0.1 -0.3 optically luminous 
(L b > 2 x 10 9 Lb.q) late-type galaxies, that 10% of the nor- 

17 Note that R mc \ has the same meaning as it did in Lehmer etal. (2007). 


mal galaxy X-ray emission originates outside of our 1 "5 ra- 
dius aperture (2 .8-6 .7 kpc at z = 0.1 -0.3) for all stacked sam- 
ples. For galaxies that were within 6f0 of more than one of the 
CDF aimpoints, we added source counts and exposure times 
from all appropriate images and exposure maps, respectively; 
however, there were very few sources in our main sample that 
met this criterion. 

Using background maps (see § 4.2 of Alexander etal. 2003 
for the ^2 Ms CDF-N and « 1 Ms CDF-S and § 4 of Lehmer 
etal. 2005b for the ^250 ks E-CDF-S) and exposure maps, 
we measured local backgrounds 7?, i oca i and exposure times 
7^. local within a 30 pixel x 30 pixel (^15" x 15") square, cen- 
tered on each source with the 1 "5 radius circle masked out. 
Here again, if a source was within 6.'0 of more than one of the 
CDF aimpoints, we summed the local backgrounds and expo- 
sure times. We estimated the expected number of background 
counts in each circular aperture 7?, by scaling the background 
counts within the square by the relative exposure times of the 
circular aperture and the square (i.e., 7?, = Tf/jocai x 7//7/.i oca l). 
This approach is similar to scaling the background counts in 


the square by the relative areas of the circular aperture and 
the square; however, by using exposure times, we are able to 
account more accurately for spatial variations in pixel sensi- 
tivity due to chip gaps, bad pixels, and vignetting. Further- 
more, comparisons between this method and the Monte Carlo 
method used in Lehmer etal. (2007) to compute £, give ex- 
cellent agreement and are most convergent for large numbers 
of Monte Carlo trials. 

When stacking galaxy populations, we excluded sources 
that were (1) classified as AGN candidates (via the crite- 
ria outlined in § 4.1), (2) within 10"of an unrelated source 
detected in the X-ray catalogs, (3) within the extent of ex- 
tended X-ray sources (see Bauer etal. 2002b, § 3.4 of Gi- 
acconi etal. 2002, and § 6 of Lehmer etal. 2005b), and (4) 
located within 3"of another late-type galaxy in our main sam- 
ple. We note that we include X-ray-detected normal galax- 
ies when stacking our samples, since we are interested in the 
average properties of the normal late-type galaxy population. 
We tested the effects of including such sources by stacking 
samples both with and without X-ray-detected sources in- 
cluded, and find similar results for both cases. Typically, the 
X-ray-detected sources contribute ^10-30% of the counts in 
the total stacked signal. For each of the subsamples of normal 
late-type galaxies outlined in Figure 4 ( gray rectangles ), we 
determined stacked source-plus-background (5 = Sf) and 
background counts {B = Bj) to determine net counts (S-B). 

For each stacked sample, we required that the signal-to-noise 
ratio [S/N = (, S-B)/y/B ; where \S-B\ and B > 20] be 
greater than or equal to 3 (i.e., ^ 99.9% confidence) for a de- 
tection. For stacked samples without significant detections, 
3cr upper limits were placed on the source counts. 

We converted the net counts obtained from each stacked 
sample to absorption-corrected fluxes and rest-frame lumi- 
nosities using a power-law SED with T = 2. Due to the fact 
that our 1 ''5 radius stacking aperture encircles only a frac- 
tion of the PSF 18 for sources at relatively large off-axis an- 
gle, we calculated aperture corrections £, for each stacked 
source i. Since we are calculating average X-ray counts from 
the summed emission of many sources of differing back- 
grounds and exposure times, we used a single, representative 
exposure- weigh ted aperture correction, £. This factor, which 
was determined for each stacked sample, was calculated as 
follows: 


. . x t , 


where T = 2/7/- The average aperture corrections (Q for 
sources in our main sample were ?^1.6, 1.8, and 1.7 for the 
SB , HB , and FB , respectively. Using our adopted SED, we es- 
timated observed mean X-ray fluxes using the following equa- 
tion: 

/ei-£ 2 =Ae x -e 1 £> ’ (6) 


where Ae x -e 2 is a bandpass-dependent count-rate-to-flux con- 
version factor that incorporates both the X-ray SED informa- 
tion as well as Galactic extinction using the column densities 
listed in § 1 . These mean X-ray fluxes were then converted 
to rest-frame luminosities following equation 4, assuming a 
photon index of T = 2. 

18 At off-axis angles 0 « 3', our If' 5 radius circular aperture contains an 
encircled-energy fraction of ^100%, 80%, and 100% for the SB, HB, and 
FB, respectively; however at 9 « 6', this fraction decreases to ^30%, 25%, 
and 25%, respecively. 



Fig. 9.— Effective photon index (r e ff) versus the logarithm of the X-ray- 
to-optical flux ratio (log/0.5-8 kev///?) for 44 stacked samples selected via 
observed properties: Lb {filled circles ), A/* (filled squares ), and SFR (filled 
triangles). The characteristic mean error bar for each quantity is given in the 
lower right-hand comer. The median logarithm of the X-ray-to-optical flux 
ratio is indicated with' a vertical dotted line (log/0.5-8 kev///? = -2.53). The 
median effective photon index for the samples that were detected in both the 
SB and HB is indicated with a horizontal dotted line (r"Jj; dian = 1.47). We 
note that all of the stacked samples with both SB and HB detections have 
r cf f > 1, consistent with normal galaxies. For the several stacked samples 
that have only SB detections, we have indicated lower-limits on T e ff. The 
shaded regions and corresponding boundaries (dashed lines) represent areas 
where X-ray-detected sources were classified as AGN candidates (for details, 
see discussion of criteria 1 and 2 in § 4.1); note that additional criteria were 
used to identify potential AGNs when generating our samples of normal late- 
type galaxies (see criteria 3-5 in § 4.1). 


5. RESULTS 
5.1. Stacking Results 

Using the stacking analysis methods discussed in § 4.2, we 
stacked the late-type galaxy samples presented in Figure 4 
(stacked samples are denoted with thick gray rectangles). 
These samples were selected using L B , M*, and SFR, which 
include 14, 17, and 13 stacked samples (44 total), respectively. 
In Table 2 {available electronically ), we tabulate our X-ray 
stacking results. We found significant (i.e., S/N £ 3) X-ray 
detections in the 0.5-2 keV and 0.5-8 keV bandpasses for all 
stacked samples. In the 2-8 keV bandpass, 15 out of the 44 
stacked samples were detected, and these samples generally 
constitute the most optically luminous and massive galaxies, 
as well as those galaxies that are most actively forming stars. 

In Figure 9, we show the effective photon index (r e ff) ver- 
sus the logarithm of the X-ray-to-optical mean flux ratio 
(log/0.5-8 keV / /r) for our 44 stacked samples {filled symbols ) 
that were selected via their observed properties. Effective 
photon indices were estimated using HB-to-SB count-rate ra- 
tios (i.e., $ 2-8 kev/$o. 5-2 kev). All 44 stacked samples have 
X-ray-to-optical flux ratios and X-ray spectral slopes consis- 
tent with normal galaxies {unshaded region in Fig. 9), sug- 
gesting that these samples are not heavily contaminated by an 
underlying population of AGNs. 19 We find a median loga- 

19 We note that two of our samples have spectral slopes that are near our 
chosen division between normal galaxies and AGNs (r e ff ~ 1). These sam- 
ples have logLfl = 10-10.5 and logA/* = 9.5-10.1 and lie at z ~ 0.66, a red- 
shift where large-scale sheets of galaxies and AGNs have been isolated pre- 
viously in the CDF-S (e.g., Gilli etal. 2003).. These galaxies have spectral 
slopes that are either statistically scattered towards T e ff = 1 (see characteris- 
tic error bar in Fig. 9) or may have some non-negligible contribution from 



FIG. 10.— Logarithm of (a) the X-ray-to-5-band mean luminosity ra- 
tio (log Lx/ L b), ( b ) the X-ray-to-stellar-mass mean ratio (log Lx/A/*), and 
(c) the X-ray-to-star-formation-rate mean ratio (logLx/SFR) versus redshift 
( filled symbols and curves) for stacked normal late-type galaxy samples se- 
lected by Lb (Fig. 4a), A/* (Fig. 4b), and SFR (Fig. 4c), respectively. For 
comparison, in each panel we have indicated the properties of the MW and 
local galaxies M101 , M82, and NGC 3256. Quoted X-ray luminosities corre- 
spond to the 0.5-8 keV bandpass and were calculated following the methods 
described in § 4.2, assuming a power-law SED with photon index of T = 2. 
Symbols and curves correspond to unique ranges of Lg , A/*, and SFR, which 
are annotated in each respective figure. For reference, in Figures 10a and 10/? 
we have plotted the corresponding values of log Lx/ L b and log Lx /A/*, re- 
spectively, for normal late-type galaxies in the local universe ( open symbols) 
using the SOI sample. In Figure 10/? we show the expected LMXB contribu- 
tion based on Gilfanov etal. (2004b; dashed line). Finally, in Figure 10c we 
show the local Lx -SFR relation and its dispersion ( dashed line with shading) 
derived by PR07 and corrected for our choice of IMF and T = 2. We note 
that roughly all of our data points in Figure 10c lie ssO.l -0.2 dex below the 
PR07 relation, which is likely due to systematic differences in how Lx and 
SFR were determined between studies. 


rithm of the X-ray-to-optical flux ratio of log / 0 . 5-8 kcv //r = 
-2.53 ( vertical dotted line), and for samples that were de- 
tected in both the HB and SB, the median effective photon 
index is T^ dian = 1.47 ( horizontal dotted line). These values 
are characteristic of galaxies dominated by X-ray binary pop- 
ulations. 

In Figure 10, we show the logarithm of the ratio of the 
0.5-8 keV luminosity (hereafter. Lx) to each physical prop- 
erty (i.e., log Lx/L b , logLx/M*, and logL x /SFR) versus red- 
shift for our samples. Each relevant quantity for the MW 
and local galaxies M101, M82, and NGC 3256 have been 
shown for comparison. We adopt 6.6 x 10 39 ergs s -1 as 
the approximate X-ray luminosity of the MW (Grimm etal. 
2002). For the local galaxies, we utilized the X-ray lumi- 
nosities from Shapley etal. (2001; hereafter SOI) for M 101 
(Lx = 1 .6 x 10 40 ergs s -1 ) and M82 (Lx = 4.6 x 10 40 ergs s -1 ) 
and Lira etal. (2002) for NGC 3256 (L x = 4.0 x 10 41 ergs s" 1 ). 
We have corrected the X-ray luminosities to be consistent 
with our use of the 0.5-8 keV bandpass and our choice of 

heavily-obscured AGNs. 


r = 2. 

For the purpose of comparing our results to those for nearby 
late-type galaxy populations, we made use of the SOI sam- 
ple of 183 normal local (D 100 Mpc) spiral and irregular 
galaxies. These galaxies were observed in the 0.2-4 keV band 
using the Einstein IPC and HRI (Fabbiano etal. 1992) and 
AGNs, with luminous X-ray emission or spectral signatures 
indicative of AGN activity, have been excised from the sam- 
ple. To avoid the inclusion of early-type SO galaxies, we chose 
to utilize 139 normal late-type galaxies from the SOI sample 
with morphological types T >2 (Sa-Irr Hubble types). 

The SO 1 sample covers ranges of L B and M* that are well- 
matched to our main sample and are representative of late- 
type galaxies in the local universe. In order to compute mean 
X-ray luminosities for samples that were directly comparable 
with our results, we divided the SOI sample into the same in- 
tervals of L fl and M* that were used for our main sample (see 
Figs. 4 a and 4b). Since several of the SOI galaxies had only 
X-ray upper limits available, we computed mean X-ray lumi- 
nosities and errors using the Kaplan-Meier estimator available 
through the Astronomy SURVival Analysis software pack- 
age (ASURV Rev. 1.2; Isobe & Feigelson 1990; LaValley 
et al. 1 992); the Kaplan-Meier estimator handles censored data 
sets appropriately. When computing these mean X-ray lumi- 
nosities, we filtered the SOI samples appropriately into dis- 
tance intervals to avoid the Malmquist bias. In Figures 10a 
and 10b, we show the corresponding values of log L x /L B and 
log Lx/M*, respectively, for the SOI sample with open sym- 
bols. By contrast, the SFRs of the local sample are generally 
too low ( <, 1-10 Mq yr -1 ) to provide a meaningful compar- 
ison with our distant 24/xm-detected sources. This is due to 
the strong positive evolution of the star-formation rate density 
with redshift (see § 1), which makes SFRs that are common 
for galaxies in our sample ( £ 1-10 M 0 yr -1 ) comparatively 
rare at z = 0. 

From Figures 10a and 10L>, it is apparent that there is signif- 
icant positive redshift evolution in lo gL x /L B and log Lx/M* 
over the redshift range of z ~ 0-1 .4. For each of the six to- 
tal selection ranges of L B and M*, the redshift progression 
of X-ray luminosities is inconsistent with a constant at the 
>99.9% confidence level. For the most optically luminous 
{L b =3-20 x 10 10 L b ,q) and massive (M* = 1-20 x 10 10 Mq) 
late-type galaxies at z = 1.4, L X /L B and L x /M* are mea- 
sured to be larger than the local values of SOI by factors of 
3.1 ± 0.7 and 9.6 ± 3. 1 , respectively. Such values are consis- 
tent with the ^(1 -Hz) 1 5-3 evolution of L£ found for the nor- 
mal late-type galaxy population, which has been constrained 
using largely the most optically luminous and massive galax- 
ies (Georgakakis etal. 2006; Ptak etal. 2007; Tzanavaris & 
Georgantopoulos 2008). 

The above results confirm the increase in L X /L B with red- 
shift found by H02. In past studies, L X /L B has been used 
as a proxy for star-formation activity (e.g., Ptak etal. 2001; 
H02; Lehmer etal. 2005), despite the fact that L B is likely 
to be somewhat sensitive to SFR (see discussion in § 3.1). 
As discussed in § 1 , L x has been shown to be strongly cor- 
related with galaxy SFR. Correlation studies of spiral galax- 
ies in the local universe have also found strong correlations 
between L x and L B such that L x oc Lg 5 (SOI; Fabbiano & 
Shapley 2002). The nonlinear relationship between X-ray 
and 5-band emission is thought to be due to the increase in 
dust obscuration with star-formation activity, which attenu- 
ates light from the 5-band more effectively than it does in 
the X-ray band. Using our sample of 47 normal late-type 



Table 3. Parametric Fitting Results For Stacked Samples 


/phys 

V 

x 2 

A 

B 

C 

0) 

(2) 

(3) 

(4) 

(5) 

(6) 

Lb 

11 

25.37 

31. 34 ±1.23 

0.85 ±0.13 

2.91 ±0.61 

AL 

14 

25.73 

33.97 ±0.87 

0.58 ±0.07 

4.48 ±0.37 

SFR 

10 

9.72 

39.91 ±0.06 

0.89 ±0.11 

0.32 ±0.88 


NOTE. — This table contains basic fitting parameters for x 2 fits to our X-ray stacking 
results. For each sample, selected by physical property / P h ys , we performed parametric 
fits for the mean 0.5-8 keV luminosity L\ following log L\ - A + B log / P h ys +C log(l +z ) . 
Col.(l): Physical parameter / p h ys used in fitting our stacking results. Col.(2): Number 
of degrees of freedom used in each fit. Col .(3): Minimum value for each fit. Col .(4)- 
(6): Best-fit values of A, B, and C with errors (90% confidence). For futher details, see 
§5.1. 


galaxies with both X-ray and 24/im detections, we found that 
Lx/ L b was positively correlated with the UV dust-extinction 
measure (Lir+LuvV^uv (Kendall’s r = 0.44; >, 99.99% con- 
fidence level), thus providing support for this hypothesis. Fur- 
thermore, it has been shown that Lx/ L b is correlated with 
^ 60 /xm/ Lioo/zm •> suggesting that more intense X-ray emission 
is associated with hotter IR colors, which is indicative of in- 
tense obscured star-formation activity. Therefore, the increase 
in Lx /L b with redshift (see Fig. 10a) observed for our galax- 
ies is likely due to an increase in their star-formation activity. 

Upon comparing Lb with M* and SFR for our galaxies, 
we find that L B is well correlated with both M* (Kendall’s 
r = 0.69) and SFR (t = 0.58); however, the scatter in the L B - 
M* relation is significantly smaller than it is for the L B - SFR 
relation (?^0.2 dex versus ^0.4 dex, respectively) implying 
that Lb traces more effectively M* rather than SFR. It is there- 
fore not surprising that we see similar redshift evolution of 
Lx/M* and Lx/ L b. 

In Figure 106, we show the estimated LMXB contribution 
to log Lx/M* ( dashed line) based on Table 5 of Gilfanov et al. 
(2004b). This value is «5-10 times lower than all mean 
values of Lx/M*, suggesting that on average LMXBs play 
a fairly small role in the X-ray emission from our stacked 
samples. Furthermore, late-type galaxies in the local uni- 
verse with similar stellar masses are often found to have 
HMXB emission that is «2-10 times more luminous than 
that expected from LMXBs (see open symbols in Fig. 106 
and Fig. 3 of Gilfanov etal. 2004c). For galaxy samples se- 
lected via SFR, we find no evidence for significant evolu- 
tion in log Lx /SFR (Fig. 10c). We note that if the majority 
of the X-ray flux from our stacked samples originated from 
X-ray-detected sources, then the AGN selection criteria 3 and 
4 from § 4.1 may potentially bias these results. However, we 
find that the more populous and unbiased X-ray-undetected 
source populations dominate the stacked 0.5-8 keV signals at 
the 80% level (see column 14 of Table 2). For each of the 
three ranges of SFR, log Lx /SFR is consistent with a constant 
value, and has a best-fit ratio of log Lx /SFR = 39.87 ± 0.07 
(X 2 = 13.1 for 12 degrees of freedom). These results sug- 
gest that the contribution from LMXBs is small and that the 
integrated X-ray emission from our late-type galaxies is dom- 
inated by HMXBs. This implies that the evolution of our late- 
type galaxy samples is likely due to changes in star-formation 
activity. 

Since the X-ray emission from our late-type galaxies is 
dominated by star-formation processes, we note that our M*- 
selected stacking results provide a relatively unobscured mea- 
sure of the star-formation activity per unit stellar mass (i.e., 


the SSFR; see Fig. 106). We find that at z « 1 the X-ray emis- 
sion per unit stellar mass is a factor of ^2-3 larger for galaxies 
with M* = 3-10 x 10 9 M© versus that observed for galaxies 
with M* = 1-20 x 10 10 M q . At z » 1 , we find that L x /M* 
is larger than its local value (SOI) by factors of 6.4 ±2.2 
and 5.8 ± 1.6 for late-type galaxies with M* = 3-10 x 10 9 Mq 
and M* = 1-20 xl0 l0 M©, respectively. These results are 
broadly consistent with observed differences in the mean SS- 
FRs found by Zheng etal. (2007) for z~ 1 galaxies of com- 
parable stellar masses, and imply that the lower-mass galax- 
ies are undergoing more significant stellar mass growth over 
z ~ 0-1 than more massive galaxies. 

In order to quantify the dependences of the X-ray luminos- 
ity on redshift and physical properties, we performed multi- 
variate parametric fitting to our stacked data. For each galaxy 
sample selected via L 5 , M*, and SFR, we investigated the red- 
shift evolution of the X-ray luminosity. For this analysis, we 
fit our data to a power-law parametric form: 

log Lx (/phys, z) = A+£log/phys + Clog(l +z), (7) 

where / p h ys is a place holder for each of the three physical 
properties (L#, M*, and SFR) used for our sample selections, 
and A, B, and C are fitting constants. For each sample, we uti- 
lized our X-ray stacking results, equation 7, and x 2 fitting to 
estimate the constants A, B, and C. For our adopted three- 
component model, we constrained A, B, and C using 90% 
confidence errors (Ax 2 = 2.7). The SOI local data points were 
not used for these fits due to differences in galaxy selection, 
instrument calibration, and AGN identification. 

In Table 3, we tabulate our constraints on x 2 > A, B, and C 
for these fits. We find that this particular choice (i.e., eqn. 7) 
of parameterization does not provide acceptable fits for galaxy 
samples selected via L B and M*. However, for galaxy samples 
selected via SFR, we find a good fit for this parameterization 
(X 2 = 9.42 for ten degrees of freedom). We constrain the evo- 
lution of logLx/SFR to be independent of or at most weakly 
dependent on redshift [pc (1 + z )0 40 ±o.53j 

Based on radio observations of distant star-forming galax- 
ies with SFR « 3-300 M© yr~ l in the CDF-N and CDF-S, 
Barger etal. (2007) reported that the X-ray upper limits for 
X-ray-undetected sources were below the level expected from 
the Lx -SFR correlaiton, thus suggesting that the correlation 
may not hold in the high-redshift universe. However, our 
stacking results suggest that the Lx -SFR correlation does 
in fact hold for average galaxies with SFR = 1-5 M© yr -1 
(SFR = 15-100 Mq yr l ) out to z « 0.5 (z « 1 .4). 

For illustrative purposes we have created Figure 1 1 , which 
shows log SFR versus log Lx for normal star-forming galax- 
ies selected from several different sources including X-ray- 
detected galaxies from our main sample {filled gray circles 
and limits ), stacked galaxies from this study {filled black 
symbols) , local galaxies from PR07 {open circles ), local ul- 
traluminous infrared galaxies from PR07 (ULIRGs; open 
squares ), and stacked z ~ 3 Lyman break galaxies (LBGs; 
stars ; see § 6 for details). For all data used in this plot, we 
have normalized SFRs appropriately to be consistent with our 
adopted Kroupa (2001) IMF and have adjusted X-ray lumi- 
nosities to correspond to the 0.5-8 keV band using a T = 2 
power-law SED. The best-fit Lx -SFR correlations for local 
galaxies from PR07 {solid curve) and z= 0-1.4 late- type 
galaxies from this study (log L x /SFR = 39.87; dotted curve) 
have been shown for reference. For comparison, we have 
shown the SFRs for the MW and local galaxies M101 , M82, 
and NGC 3256. 



log L x (ergs s ‘) 

FlG. 11. — Logarithm of the star-formation rate log SFR versus the loga- 
rithm of the X-ray luminosity log Lx for normal late-type galaxies. X-ray- 
detected sources from our main sample are indicated as small gray dots. 
Sources that were not detected in the 0.5-8 keV bandpass but were detected 
in either the 0.5-2 keV or 2-8 keV bandpasses are shown as upper limits. 
Results from our X-ray stacking analyses of late-type galaxies selected via 
observed SFR are shown as large filled circles, squares, and triangles, which 
have the same meaning as in Figure 10c; z ~ 3 LBGs that were both un- 
corrected and corrected for AGN contamination have been shown plotted as 
a filled star and open star, respectively. For comparison, we show the local 
galaxy sample from PR07, which includes normal late-type galaxies (open 
circles) and ultraluminous infrared galaxies (ULIRGs; open squares)-, the 
best-fit PR07 relation is shown as a solid curve. The SFRs for the MW and 
local galaxies Ml 01 , M82, and NGC 3256 are indicated (horizontal dashed 
lines ) . 



Fig. 12.— Cumulative X-ray luminosity dependent AGN fraction 
£l(L 2 - 8 kev) versus L 2 - 8 keV for late-type galaxies with Lb ~ 2 x 10 10 Lb.q 
at z £ 0.6 (filled circles) and z ~ 0.6-1 (open triangles). For reference, we 
have indicated the corresponding AGN fraction for Lyman break galaxies at 
z ~ 3 (filled stars-, see § 5 for further details). We have indicated the median 
X-ray detection limit for galaxies in each redshift range (downward-pointing 
arrows along the x -axis). The solid curve represents the best-fit relation for 
#l(L 2 _h k C v), which was fit using all late-type galaxies with z ss 0-1 and 
Lb 2 x 10 10 Lfi.Q in our main sample. For reference, we have shown the 
estimated AGN fraction for galaxies with Lb ~ 5 x 10 9 , 10 10 , and 10 1 1 Lg, © 
(dotted curves: see § 5.2 for further details). 



Redshift ( z ) 

Fig. 13.— Cumulative redshift-dependent AGN fraction g z (z) ver- 
sus redshift for late-type galaxy samples with Lb £ 10 10 Lb,q and 
L 2 _g kev ~ 10 41 5 ergs s" 1 . We find no significant evolution of g z (z) over 
the redshift range z ~ 0.1 -0.8. 


5.2. AGN Contribution to Stacked Signals 

In this section we estimate the contribution to our stacked 
signals from contaminating AGNs that have X-ray luminosi- 
ties below our X-ray detection limit. This analysis is sim- 
ilar in nature to that in §§ 3.1 and 3.2.2 of Lehmer etal. 
(2007), which was performed for early -type galaxies. We im- 
plement the observed cumulative AGN fraction fc : the frac- 
tion of galaxies harboring an AGN with 2-8 keV luminosity 
of L 2-8 kev or greater. Hereafter, we compute fc by taking 
the number of candidate AGNs in a particular galaxy sample 
with a 2-8 keV luminosity of L 2 ~ % k e v or greater and dividing 
it by the number of galaxies in which we could have detected 
an AGN with luminosity L 2 ~ 8 kev • The latter number is com- 
puted by considering the redshift of each galaxy and its corre- 
sponding sensitivity limit, as obtained from spatially varying 
sensitivity maps (see § 4.2 of Alexander et al. 2003 and § 4 of 
Lehmer etal. 2005b); these sensitivity maps were calibrated 
empirically using sources detected by wavdetect at a false- 
positive probability threshold of 1 x 10“ 5 . 

The quantity fc is not only a function of L 2 -% kev > but is also 
dependent on the selection of the galaxy sample considered: 
in our case redshift and the physical property / P h ys (i.e., L B , 
M* , and SFR) of the galaxy population may plausibly play a 
role in f c . In our analyses, we assume that each of the respec- 
tive dependencies (i.e., L 2 -8 kev, z, and / p h ys ) are independent 
of each other, such that f c oc gdU-% kev) x g z {z) x g p {f phy s ), 
where g L , g-, and g p represents the functional dependence 
of the cumulative AGN fraction for each indicated variable 
£ 2-8 kev, z , and / P h ys , respectively. We made use of the 
2-8 keV bandpass because of its ability to probe relatively 
unattenuated X-ray emission in a regime of the X-ray spec- 
trum where we expect there to be minimal emission from nor- 
mal galaxies (see also criterion 1 of § 4.1 for further details). 
In total 62 (^51%) of our 121 X-ray-detected AGN candi- 
dates had 2-8 keV detections; we use these AGNs in our AGN 
fraction analyses. 

We began constructing fc by estimating the shape 
of g L (L 2-8 kev) using late-type galaxies with z = 0-1 and 
L B ^ 2 x 10 10 L b q , an optical luminosity regime where we 
have a relatively large number of sources and are sufficiently 
complete out to z = 1 (see Fig 4a). We split this sam- 






Fig .14.— Cumulative AGN fraction g p {f p hys ) versus (a) Lb , ( b ) M* , and (c) SFR for late-type galaxies with L 2-8 keV £ 10 41 5 ergs s 1 at z < 0.6 {dashed curves) 
and z = 0.6-1 ( dotted curves). In each plot, we have indicated the best-fit relation {solid line), which is calculated using all galaxies with L 2-8 keV £ 10 415 ergs s -1 
and z < 0.6. The downward-pointing arrows indicate which values of / p h ys were chosen for computing g 0 in equation 8. For reference, the bin sizes of / p h ys used 
to compute each value of g p {f p by&) and the typical errors of g p {f p hy&) have been indicated in the lower right-hand comer of each plot. 


pie into two subsets about z = 0.6, to test whether there 
is substantial evolution in the shape and normalization of 
gifL 2-8 kev) over this redshift range. In Figure 12, we 
show gi(L 2 _ 8 kev) for galaxies in the redshift ranges z £ 0.6 
(^median =0.51; filled circles) and z = 0.6-1 (z me dian =0.84; 
open triangles). From Figure 12, we see that the overall shape 
and normalization of gi{L 2 _ 8 kev) for late-type galaxies with 
L B > 2x 10 10 L 5)0 is similar for galaxies at Zmedian = 0.51 
and Zmedian = 0.84. We fit the shape of g L (L 2 _ 8 ke v) using 
least-squares fitting of the 2-8 keV luminosity dependent cu- 
mulative AGN fraction using all galaxies from z = 0-1 with 
L B Z 2 x 10 10 L b .q. For these fits, we found that the data 
were well-fit by an exponential function, which we parameter- 
ized as logg L (L 2 _ 8 kev) oc aexp[-6(logL 2 _ 8 kev - 39) 2 ] , where 
a and 6 are fitting constants. By construction, this function is 
only valid for logL 2 _s kev > 39, which is « 1-2 orders of mag- 
nitude less luminous than a typical stacked X-ray luminosity 
of our late-type galaxy samples (see Table 2 and Fig. 10). 
We find best-fit values of a = -0.8 and b = -0.05; in Fig- 
ure 12 ( thick black curve), we show our best-fit relation for 
gL(^2-8 kev)- 

We constrained further the redshift evolution of fc [i.e., the 
shape of g z (z)] by dividing our main late-type galaxy sample 
into five nearly independent redshift bins (with z = 0.1 -0.8) 
and calculating fc for fixed ranges of L 2 _ g k C v and / P h ys . In 
Figure 13, we show g-(z) as a function of redshift for late-type 
galaxies with 8 ke v £ 10 41 5 ergs s" 1 and L B £ 10 10 Lj* s0 , 
the approximate completeness limit at z « 0.8 (see Fig. 4a). 
Using these data and \ 2 fitting [assuming a ( 1 + z) n depen- 
dence], we constrained the redshift evolution of g z (z) to be 
proportional to (1 +z) a,5±a97 ; similar results were found for 
different ranges of L 2 _ 8 kev and Lb . This result differs from 
that observed for early-type galaxies, where the AGN fraction 
and mean AGN emission has been found to evolve as «( 1 +z) 3 
(e.g., Brand etal. 2005; Lehmer etal. 2007). 20 We therefore 
conclude that there is little redshift evolution in the late-type 

20 We note that the lack of evolution of the AGN fraction does not nec- 


galaxy AGN fraction over the redshift range z = 0.1 -0.8, and 
hereafter we assume that g z (z) remains roughly constant out 
to z= 1.4. 

To constrain the overall dependence of fc on / P h ys [i.e., 
g/,(/ p h y s)], we calculated the cumulative AGN fractions for 
late-type galaxy samples with / p h ys by holding the ranges 
of La - 8 kev and z fixed and varying / p h ys . In Figure 14, we 
show g p (/ph ys ) versus logL 5 (Fig. 14a), logM* (Fig. 146), 
and logSFR (Fig. 14c) for L 2 _ 8 kev ~ 10 41 5 ergs s" 1 at z £ 0.6 
(1 dashed curves) and z = 0.6-1 ( dotted curves). We calculated 
g P (f phys) for intervals of / P h ys where we are approximately 
complete at z = 0.6 (for the z ^ 0.6 interval) and z = 1 (for the 
z = 0.6-1 interval; see Fig. 4). For each sample, we again uti- 
lized least-squares fitting to approximate the / P h ys dependence 
of g P (/ P h ys ). These fits were performed using all data over the 
redshift range z = 0-1 assuming a functional dependence of 
l°gg/>(/phys) oc clog/phys. In each panel of Figure 14, we show 
the best-fit solutions for g p (/ p h ys ) with the gray lines. We find 
that the AGN fraction is strongly dependent on L B (Fig. 14a) 
and M* (Fig. 146), such that more optically-luminous and 
massive galaxies have larger AGN fractions; this result is 
consistent with other studies of the AGN host galaxies (e.g., 
Kauffmann etal. 2003b; Nandra etal. 2007; Silverman etal. 
2007). Also, the AGN fraction seems to be mildly dependent 
on the galaxy SFR; however, this is likely due to the fact that 
the SFR is larger on average for more massive galaxies. 

Based on the above estimates of the shapes of gL^i-% kev), 
g z (z), and g p (f phys), we approximated empirically f c for each 
choice of / phys following: 

log fc = log go+a exp[-b(logL 2 -8 kcV-39) 2 ] +clog/ phys , (8) 
where go represents the normalization of each relation based 

essarily imply that the optically-luminous (i.e.. Lb £ 10 10 Lb.q) late-type 
galaxy AGN number density is constant over z ~ 0-1 . For example, LJ has 
been shown to fade by ss 1 mag since - a 1 (e.g., Lilly etal. 1995; Wolf etal. 
2003; Faber etal. 2007), suggesting that there are fewer optically luminous 
late-type galaxies in the local universe than at z = 1 . 



upon the value of fc at £ 2-8 kev > 10 41 5 ergs s 1 and 
L b = 10 10 £*,©, M* = 10 10 A#©, and SFR = 2 M© yr' 1 , for 
galaxy samples selected via £#, M*, and SFR, respectively. 
We note that only go and c are dependent upon our choice 
of /physi we find for the set of physical parameters / p h ys = [L B , 
M*, SFR] that g 0 = [6. 1 x 1(T 14 , 8.2 x 10" 8 , 0.20] and c = [ 1 .2, 
0.67,0.66]. For reference, in Figure 12 we have shown curves 
of fc for L b = 5 x 10 9 , 1 0 10 , and 10 11 £#.© (< dotted curves). 

To estimate the AGN contamination expected for each of 
our stacked samples presented in § 5.1, we followed closely 
the procedure in § 3.2.2 of Lehmer etal. (2007). For com- 
pleteness, we outline this procedure below. 

For each stacked sample, we used equation 8 to compute 
cumulative AGN fractions fc. In Figure 15 a, we show fc 
for two of our £ fl -selected samples: z ^ 0.2 galaxies with 
L b = 10 95 -10 10 ( 1 dashed curve) and z ~ 1.3 galaxies 

with L b = 10 10 5 — 10 1 1 3 £#© ( solid curve). We then converted 
fc for each sample to a differential AGN fraction fo (i.e., 
the 2-8 keV luminosity dependent fraction of galaxies har- 
boring AGNs within discrete X-ray luminosity bins of width 
A log £ 2-8 kev = 0.5; see Fig. 15b). For each sample, we cal- 
culated the 2-8 keV luminosity dependent fraction of galaxies 
that were below our £ 2 -8 kev detection limit f B (Fig. 15c). We 
then multiplied f D and f B to estimate the 2-8 keV luminosity 
dependent fraction of galaxies harboring an AGN that was un- 
detected in the Chandra exposures fc (i.e., /(/=/ D x f B \ see 
Fig. 15 d)\ these AGNs would not have been removed for our 
stacking analyses. Finally, we approximated the total 2-8 keV 
AGN contamination £ 0 . 5-2 kev(contam) of each stacked sam- 
ple using the following summation: 

£2-8 kev(contam) = X^ 2 - 8 keV,/, ( 9 ) 

i 

where the summation is over all bins of A log £ 2-8 kev = 0.5 
in the range log£ 2 -8 kev = 39-45. We converted each value of 
£ 2-8 kev(contam) to estimated values of £ 0 . 5-8 kcv(contam) and 
£0 5-2 kev(contam) by assuming the contaminating AGN emis- 
sion roughly follows an X-ray SED described by a power-law. 
In order to constrain the average photon index of the power- 
law, we stacked all X-ray-detected AGNs in our stacked sam- 
ples with 0.5-8 keV luminosities below 10 42 ergs s -1 . For 
these AGNs, we find a stacked effective photon index of 
Teff = 0.97 ± 0.03. If the intrinsic value of the photon index is 
T = 2, at the median redshift of our main sample (z me dian = 0.5), 
r e ff ~ 1 corresponds to an intrinsic X-ray column density 
of Ah ~ 1-2 x 10 22 cm -2 . If we assume r e ff = 1 describes 
well the effective SED of the X-ray-undetected AGNs in our 
stacked sample, we find that AGN contamination can account 
for ^1-30% (median of &5%) of the 0.5-8 keV emission 
from our stacked samples, suggesting that AGNs are not pro- 
viding a significant contribution to our stacked results. We 
note that the X-ray SED used in this calculation has an impor- 
tant effect on the overall estimate of the AGN contamination. 
Since our estimates for contamination in the 0.5-8 keV and 
0.5-2 keV bandpasses decrease as T e ff decreases, the amount 
of contamination in our samples may be underestimated if our 
choice of T c ff is too flat; however, we find that for conserva- 
tive choices of T c ff (i.e., r e ff 2) that are representative of 
even unobscured AGNs (for reference, see Fig. 5), the AGN 
contamination remains low ( <, 40%) and has no material ef- 
fect on our results. 

It is also worth noting that if there exists a large popula- 
tion of intrinsically luminous yet heavily obscured AGNs in 
our high-redshift galaxy samples that are X-ray-undetected in 



l°g £ 2 -b kev 

Fig. 15.— (a) Cumulative AGN fraction (i.e., the fraction of 
galaxies harboring an AGN with a 2-8 keV luminosity of £ 2 -g kev 
or greater), /c, versus log £ 2-8 kev for two of our stacked samples: 
Z ~ 0.2 galaxies with L B = 10 95 -10 10 £*,© {dashed curve) and 
1.3 galaxies with L B = 10 10,5 -10 11,3 £*,© {solid curve). Each 
curve was computed following equation 8. (b) Differential AGN 
fractions (i.e., the fraction of galaxies harboring an AGN in discrete 
bins of width Alog£ 2 -g kev = 0.5), f D , versus log£ 2 _ 8 kev- ( c ) Frac- 
tion of late-type galaxies for which we could not have detected an 
AGN with a 2-8 keV luminosity of £ 2 - 8 kev, /s, versus log£ 2 _ 8 ke v. 
{d) Fraction of galaxies harboring AGNs in our optically luminous 
faded samples that would remain undetected due to sensitivity lim- 
itations, fc = fo x f B , versus log£ 2 _g k c v; these galaxies would not 
have been removed from our stacking analyses. 


the CDFs, then our conclusions above could be somewhat dif- 
ferent. If such a population were present and had significant 
influence on our stacked results, then we would expect to find 
stacked X-ray spectra that were relatively flat (r e ff £1); how- 
ever, as discussed in § 5.1 , we find that very few of our stacked 
samples could have such flat spectra in our stacked samples, 
suggesting that such a population, if present, does not have a 
strong effect on our results. 

In Table 2 (cols. 25-27), we have provided the estimated 
fractional AGN contribution to each stacked sample [i.e., 
£Ei-£ 2 (contam)/££ 1 _£ 2 ] for the FB, SB, and HB using the 
technique described above and an assumed T C ff = 1 . We find 



that the estimated AGN contamination is most significant for 
galaxy samples with large values of L B , M*, and SFR. 

6. EXTENSION TO DISTANT LYMAN BREAK GALAXIES 

As noted in § 1 , the global star-formation rate density has 
been observed to increase with redshift out to z ~ 1-1 .5. In- 
vestigations of the most-distant LBGs at z ~ 2-7 show that 
the star-formation density peaks around z ~ 1-3 and gradu- 
ally declines toward higher redshifts (e.g., Steidel et al. 1999; 
Giavalisco 2002, 2004b; Bouwens etal. 2004, 2005; Dickin- 
son etal. 2004). To investigate whether the mean X-ray ac- 
tivity from normal late-type galaxies follows a similar trend, 
we study the X-ray properties of a sample of z = 3.01 ± 0.24 
LBGs, which were identified as U-band “dropouts” through 
the GOODS project (see Lehmer et al. 2005a for details). 

We filtered the GOODS z ~ 3 LBG sample (Lee et al. 2006) 
to include only LBGs that (1) were within the central ^4f0 of 
the « 2 Ms CDF-N and ^1 Ms CDF-S, and (2) had rest-frame 
B-band luminosities that were similar to the most optically lu- 
minous ( L b = 3-20 x 10 10 L B ©) z = 0-1 .4 late-type galaxies 
used in this study. B-band luminosities were calculated by 
applying ^-corrections to the zsso-band flux (from GOODS), 
where /c-corrections were derived using an SED appropriate 
for LBGs (see § 2.2 of Lehmer etal. 2005a for details). We 
found that 85 z~ 3 LBGs from the Lehmer et al. (2005a) sam- 
ple satisfied these two selection criteria. 

We identified three X-ray-detected LBGs, which all had 
L 2-8 kev £ 10 43 ergs s’ 1 . Due to their high X-ray luminosi- 
ties, we classified these sources as obvious AGNs. After 
removing these three X-ray-detected AGNs from our z ~ 3 
LBG sample, we performed X-ray stacking analyses as de- 
scribed in § 4.2. We found a significant (3.1a) detection in 
the 0.5-2 keV bandpass, which corresponds roughly to rest- 
frame 2-8 keV emission. Assuming an intrinsic power-law 
X-ray spectrum with a photon index of T = 2, we found a mean 
0.5-8 keV luminosity of L\ = (4.9 ± 1 .4) x 10 41 ergs s _1 for 
our z ~ 3 LBGs, a value that agrees well with previous stud- 
ies (e.g., Lehmer etal. 2005a; Laird etal. 2006). 

Since our LBGs reside in the high-redshift universe, AGN 
contamination is expected to have a more significant ef- 
fect on these results than it did for our z = 0-1 .4 late- type 
galaxies. The median X-ray luminosity detection limit is 
~8 x 10 41 ergs s _1 and ^10 43 ergs s _1 for the 0.5-2 keV and 
2-8 keV bandpasses, respectively. To estimate the AGN con- 
tamination, we followed the approach outlined in § 5 .2, which 
made use of the 2-8 keV luminosity dependent AGN fraction 
(/c). In Figure 12 (filled stars), we show the X-ray luminos- 
ity dependent AGN fraction £llbg(£ 2-8 kev) for z ~ 3 LBGs. 
We note that at L 2 ~ 8 k e v ~ 10 43 ergs s _1 , the cumulative AGN 
fraction for z ~ 3 LBGs is a factor of ^3-4 times larger than 
that computed for our z = 0-l late- type galaxy sample with 
similar optical luminosities. Using the functional form for f c 
presented in equation 8, but with a ^3-4 times larger normal- 
ization factor for z ~ 3 LBGs (i.e., go.LBGs ~ 3-4 go), we find 
that AGN emission may plausibly account for «50-70% of 
the stacked 0.5-2 keV counts. 

In Figure 16, we show the X-ray-to-B-band mean lumi- 
nosity ratio for L B = 3-20 x 10 10 L b .q star-forming galaxies 
(i.e., late-type galaxies and LBGs) as a function of the age of 
the Universe. Together, these data span ?^85% of cosmic his- 
tory (i.e., out to z ~ 3). As presented in § 5.1, L\/L B shows 
significant evolution over the redshift range z = 0-1 .4, and af- 
ter correcting for AGN contamination,, we find that Lx/L B is 
similar for z ~ 3 LBGs and z = 1 .4 late-type galaxies. This 
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Fig. 16.— Logarithm of the X-ray-to-fi-band mean luminosity ratio 
log L\ /Lb versus the age of the Universe (for our adopted cosmology, the 
current age of the Universe is 13.47 Gyr) for star-forming galaxies with 
Lb = 3-20 xlO 10 Lb.q . For reference, redshift has been plotted along the 
top axis. We have included mean values of L\/L B for the SOI local sam- 
ple of late-type galaxies ( open circles), our stacked samples at z ~ 0.5-1 .4 
(large filled circles ), and z ~ 3 LBGs that were both uncorrected (filled star) 
and corrected (open star) for AGN contamination. At z « 0.1 -0.4, we have 
plotted Lx /L b for individual late-type galaxies from our main sample (small 
filled circles and upper limits); the mean values and errors for these galaxies, 
computed using ASURV, has been indicated. 


result suggests that the non-AGN X-ray emission for the most 
luminous star-forming galaxies may flatten near z~ 1.4-3, 
which has been predicted roughly from simulations of how the 
normal-galaxy X-ray emission is expected to respond due to 
global changes in the star- formation rate density (e.g., Ghosh 
& White 2001). 

To test whether the Lx-SFR correlation is similar for 
z ~ 3 LBGs as we found for late-type galaxies at z = 0-1 .4, 
we approximated absorption-corrected SFRs for the z ~ 3 
LBGs using UV band emission. These SFRs were ap- 
proximated following SFR = 9.8 x 10~ n 7Luv, where 7^6 
is the absorption-correction factor (see Giavalisco etal. 
2004b). As described in § 3.3, we approximated the UV 
luminosity following L uv = 33^/^(2800 A); however, here 
/*,(2800 A) was derived using our adopted LBG SED. We 
find that the mean SFR for optically-luminous z ~ 3 LBGs 
is ?s60 M© yr _1 . After correcting the mean stacked X-ray 
luminosity for AGN contamination, we find an X-ray-to- 
SFR ratio of logLx/SFR « 39.4-39.6, which is suggestively 
lower than its value for z = 0-1.4 late-type galaxies (i.e., 
logLx/SFR = 39.87 ±0.04). We note, however, that the esti- 
mation of the mean SFR for these galaxies is based solely on 
the 2800 A emission and is therefore highly uncertain; further 
constraints on the mean infrared luminosity of these sources 
would help considerably. For reference, in Figure 1 1 we have 
plotted log SFR versus log Lx for the z ~ 3 LBGs both uncor- 
rected (filled star) and corrected (open star) for AGN contam- 
ination. 

It is interesting to note that the AGN-corrected LBG sam- 
ple appears to have Lx and SFR values similar to those of 
local ULIRGs, which are relatively X-ray underluminous for 
their derived SFRs. PR07 suggest that galaxies with SFRs that 
begin to exceed £ 50-100 M© yr -1 , similar to the ULIRGs, 
may have X-ray properties that are completely dominated by 
HMXBs. In such systems, other X-ray-emitting populations 
that may normally be significant in more quiescent galax- 



ies (e.g., LMXBs, hot gas, supemovae and their remnants, 
etc.) would collectively contribute only a negligible fraction 
of the total X-ray emission, thus causing the overall galax- 
ies to appear X-ray underluminous compared with the L\- 
SFR correlation. Another possibility is that as the absorp- 
tion within the most actively star-forming galaxies increases 
with SFR (see discussion in § 5.1), the X-ray emission from 
point-source populations may become significantly obscured. 
Deep X-ray observations of local LIRGs and ULIRGs with 
Chandra and XMM-Newton have found significant X-ray ab- 
sorption in their point-source populations (see, e.g., Lira et al. 
2002; Zezas etal. 2002, 2003) thus providing some support 
for this possibility. Future investigations into the nature of the 
X-ray populations of ULIRGs could help resolve this issue. 

7. SUMMARY AND FUTURE WORK 

We have investigated the X-ray emission from 2568 nor- 
mal late-type galaxies over the redshift range z = 0-1 .4 that lie 
within the Chandra Deep Fields (CDFs). Our late-type galaxy 
sample was constructed primarily using color-magnitude dia- 
grams, which incorporated rest- frame U -V color and abso- 
lute V-band magnitudes, to isolate blue late-type galaxies (see 
§ 2 for details). In total, 225 (^9%) of our late-type galax- 
ies were detected individually in the X-ray band. Based on 
X-ray and optical spectral properties, X-ray-to-optical flux 
ratios, the correlation between X-ray luminosity and star- 
formation rate, and comparisons between infrared and radio 
properties, we infer that 121 («53%) of the X-ray-detected 
late-type galaxies are dominated by AGN emission. The re- 
maining 104 X-ray-detected galaxies had X-ray and multi- 
wavelength properties consistent with those of normal late- 
type galaxies with X-ray emission dominated by X-ray bina- 
ries (HMXBs and LMXBs). To study the X-ray emission and 
evolution from large representative populations of late-type 
galaxies (i.e., including the X-ray-undetected sources), we 
utilized X-ray stacking analyses of galaxy populations with 
AGN candidates removed. We stacked normal-galaxy sam- 
ples that were selected via their rest-frame 5-band luminosity 
(Lb), stellar mass (M*), and star- formation rate (SFR) in red- 
shift bins (see § 5.1). Furthermore, we compared these results 
with those found for a sample of z ~ 3 LBGs from Lehmer 
etal. (2005a). In the points below, we summarize our key 
findings: 

1 . We obtained significant detections in the 0.5-2 keV and 
0.5-8 keV bandpasses for all of our stacked samples. We es- 
timated that LMXBs and low-level AGNs provide only low- 
level contributions to the stacked X-ray emission from our 
samples and that HMXBs constitute the dominant X-ray- 
emitting component. Therefore, for these galaxies the X-ray 
emission is tracing primarily star-formation activity. Nor- 
mal late-type galaxy samples selected via Lb and M* show 
significant (at the >99.9% confidence level) evolution in 
their average X-ray properties from z = 0 to 1.4. For the 
most optically luminous (Lb ~ 5 x 10 10 Lb.q) and massive 
(M* « 3 x 10 10 M©) late-type galaxies at z = 1.4, Lx/ L b and 
Lx/M* are measured to be larger than their local values by 
factors of ^2-4 and «7-13, respectively. 

2. We find that late-type galaxies of lower stellar mass gen- 
erally have larger X-ray-to-stellar-mass mean ratios (Lx/M*) 
than their higher-mass analogs. Over z « 0.2-1, galaxies 
with M* « 6 x 10 9 M 0 are a factor of ^2-3 times more 
X-ray luminous per unit stellar mass than galaxies with M* ~ 
3 x 10 10 M o . 


3. We characterized the X-ray properties of 888 24/im- 
detected late-type galaxies in our sample. The 24fim data 
allowed us to select galaxy samples selected via their SFRs. 
For these samples, we found that the X-ray luminosity is well- 
predicted by a constant Lx-to-SFR ratio, similar to the L\- 
SFR correlation reported by previous authors (e.g., PR07). 
This implies that the Lx-SFR correlation holds out to at least 
z = 0.5, 1 , and 1 .4 for galaxies with SFR « 2, 10, and 50 M© 
yr -1 , respectively, and supports the idea that the strong X-ray 
evolution observed for normal late-type galaxies selected via 
Lb and M* is likely due to strong changes in SFR. 

4. The X-ray properties of our most optically-luminous 
(L b = 3-20 xl(r° L b ,q) late-type galaxies at z = 1.4 are 
comparable to those for z ~ 3 LBGs with similar optical 
luminosities, once X-ray-undetected AGN contamination in 
the LBG population has been accounted for. This suggests 
that there may plausibly be a flattening in the Lx/L^-z 
relation for optically-luminous star-forming galaxies between 
z ~ 1 .4-3 . We estimate a mean SFR of ?«60 M© yr -1 for these 
LBGs. We find that the observed mean X-ray luminosity is 
suggestively underluminous based on the Lx-SFR correlation 
prediction; this result is similar to that found for local 
ULIRGs with comparable SFRs. 

The above results can be improved greatly through (1) the 
study of late-type galaxy populations in other existing multi- 
wavelength extragalactic survey fields that contain Chandra 
observations, (2) additional observations of the CDFs to even 
better sensitivity levels than are currently available, and (3) 
observations with future X-ray missions with imaging capa- 
bilities that are complementary to those of Chandra. 

The areal footprint of the CDF regions used in this study 
totals ^0.18 deg 1 2 (see § 2.1 for details), which severely lim- 
its the number of galaxies being studied at z £ 0.2. The red- 
shift range z = 0-0.2 spans «2.5 Gyr of cosmic look-back 
time, compared with the ^6.6 Gyr over z = 0.2-1 .4, that we 
are most effectively studying above. Several complemen- 
tary deep and wide extragalactic Chandra survey fields that 
have recently been conducted or are in progress can improve 
the present situation. Three such ideal survey fields are the 
^200 ks All- wavelength Extended Groth Strip International 
Survey (AEGIS; e.g., Nandra etal. 2005; Davis etal. 2006), 
the «50 ks Chandra Cosmic Evolution Survey (C-COSMOS; 
PI: M. Elvis; see also, Scoville etal. 2007) , 21 and the ^5 ks 
NOAO Deep Wide-Field Survey (NDWFS; e.g., Murray et al. 
2005). 

Our study shows that identifying and removing contaminat- 
ing AGNs is a crucial ingredient to studying the stacked X-ray 
properties of normal-galaxy samples (see § 5.2). In order to 
avoid significant levels of AGN contamination in shallower 
wide-area Chandra surveys such as AEGIS, C-COSMOS, and 
NDWFS, X-ray studies of normal galaxy populations must 
be limited to lower redshift intervals. If we require a de- 
tection limit of Lq. 5-2 keV = 10 41 5 ergs s -1 over regions of 
each survey where the Chandra sensitivity is optimal (i.e., 

21 We note that the C-COSMOS observing strategy has been designed to 
have several overlapping «50 ks ACIS-I exposures; in these regions, the total 
exposure will reach «200 ks due to the overlaping exposures. However, the 
analyses presented in this paper utilize only high-quality Chandra imaging at 
off-axis angles of /S 6'0 (see § 2.1 for additional justification), where there 
is very little overlap in the C-COSMOS Chandra pointings. We therefore 
do not consider the overlapping regions of the C-COSMOS exposures in this 
discussion. For additional information regarding C-COSMOS, see http://cfa- 
www.harvard.edu/hea/cos/C-COSMOS .html. 


<, 6.'0 from the Chandra aim-points), then we estimate that 
AEGIS, C-COSMOS, and NDWFS could effectively be used 
to study normal late-type galaxy populations at redshifts less 
than ~0.5, «0.25, and «0.1 , respectively. When factoring in 
the areal coverage of these fields, we estimate that AEGIS and 
C-COSMOS will contain ^250-350 late-type galaxies with 
L b ~ 3 x 10 9 L b © at z - 0. 1-0.2 (i.e., a factor of larger 
than the CDFs) and that AEGIS, C-COSMOS, and NDWFS 
taken together will have « 100-250 similar galaxies at z £ 0. 1 
(i.e., a factor of ^22 times larger than the CDFs). Study- 
ing these galaxy populations would significantly improve con- 
straints on the X-ray properties of z ~ 0.2 late-type galaxies. 

Deeper Chandra observations of the CDFs would provide 
additional insight into the X-ray properties of the normal 
late-type galaxy populations presented in this study and en- 
able us to extend our analyses to higher redshifts. Typi- 
cal galaxies in our sample have mean X-ray luminosities of 
Lx « 10 39 5 — 10 40 ergs s -1 . In the ^2 Ms CDF-N, where our 
Chandra sensitivity is greatest, we expect that galaxies with 
Lx « 10 39 5 ergs s -1 should be detectable out to z ~ 0.2. Out 
of the 1 1 late-type galaxies with L B ^ 10 9 Lfl ; © at z ^ 0.2, we 
find X-ray detections for 8 of them (?^72%). Deeper Chan- 
dra observations over the CDFs will allow for the detection 
of such galaxies out to progressively higher redshifts.. For 
Chandra exposures of «5 Ms and ^10 Ms, a source with 
Lx « 10 39 5 ergs s _1 should be detectable out to z ~ 0.3 and 
z ~ 0.4, respectively. If we pessimistically assume an X-ray 
detection fraction of ^70%, we estimate that ^30 and ?^50 
normal galaxies would be detected individually for each re- 
spective deep exposure. The individual detection of these 
galaxies would allow for improved constraints on both the 
X-ray luminosity functions of normal late-type galaxies at 
higher redshifts and the low-luminosity AGN contributions 
to our stacked signals. At higher redshifts, an ^5 Ms and 
^10 Ms Chandra exposure would enable us to effectively 
perform studies of normal late-type galaxies, similar to those 

22 For additional information regarding the future X-ray mis- 
sions XEUS and Generation-X , see http://www.rssd.esa.int/XEUS/ and 


presented in this paper, out to z ~ 2 and z « 2.5, respectively. 
Such a data set would provide, for the first time, a reliable 
constraint on the X-ray emission from normal star-forming 
galaxies near the peak of the global star-formation rate den- 
sity at z ~ 1 .5-3 . 

In addition to the improvement that additional Chandra 
observations could provide, future X-ray missions such as 
XEUS and Generation-X 22 should enable new scientific in- 
vestigations of distant normal galaxies. XEUS will be able 
to place significant spectral constraints for sources with 
0.5-2 keV fluxes down to ^ 10“ 17 ergs cm -2 s -1 , a level fainter 
than many of the average fluxes derived from our samples. 
Generation-X is planned to provide imaging with ^0''1 per 
resolution element and will easily probe to flux levels of 
«10~ 19 ergs cm -2 s -1 (Lx ~ 5 x 10 3 * ergs s -1 at z = 1). At 
these levels, the first detailed investigations of the evolution 
of the normal late-type galaxy X-ray luminosity function can 
be performed effectively out to z £ 1 without confusion prob- 
lems due to the crowding of large numbers of sources. Fur- 
thermore, these observations will allow for new constraints to 
be placed on the X-ray populations of all normal galaxies in 
the observable universe that are offset by more than ?^0.8 kpc 
from their host-galaxy nuclei. 
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